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We study physics on p(770) and u>(782) in framework of chiral constituent quark model. The 
effective action is derived by proper vertex method, which can capture all order information of chiral 
expansion. The N^ 1 expansion is also studied systematically. It is shown that the momentum 
expansion at vector meson energy scale converges slowly, and the loop effects of pseudoscalar meson 
play an important role at this energy scale. We provide a method to prove the unitarity of S-matrix 
■ in any low-energy effective theory of QCD. Phenomenologically, we study decays for p — > nn, 

p — > e + e~ , oj — > 7r + 7r~, p — > 771-* and uj —* p° — u mixing and their mass splitting, pion 

form factor, 1 = 1 = 1 phase shift and light quark masses at vector meson energy scale. These 
results include all order contribution of vector meson momentum expansion and is up to next to 
leading order of iV" 1 expansion. All of these theoretical predictions agree with data very well. 
The unitarity of S-matrix yielded by this framework is examined. The Breit-Winger formula for 
resonance propagator is derived. 
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O ■ !• INTRODUCTION 

O , 

In past thirty years, the effective theory studies on physics of light vector meson resonances(p, u>, K* and <p) have 
Oh| attracted much interests. Works in this field have helped us to understand various aspects on low energy QCD and 
properties of light vector mesons. A challenging subject in these studies is to construct a well-defined chiral expansion 
theory at vector meson scale. In ref. jl| authors have developed a chiral expansion theory including vector mesons. 
In this framework one can successfully deal with two- vector-meson processes (one vector meson is in initial state and 
one is in final state) However, it was also pointed out in ref. Q that this framework is failed to study single- 
vector-meson processes. So far, the most studies on single- vector-meson processes are still limited at the lowest order 
of vector meson four- momentum expansion and N^ 1 expansion. Theoretical reason is that convergence of momentum 
expansion in these processes is unclear, and technical reason is that the calculation beyond the leading order is very 
difficult. For example, method of chiral perturbative theory(ChPT) j| is impractical for these processes because more 
and more free parameters are required as perturbative order raising. Obviously, this is a serious problem that has 
to be resolved, since it makes the models' calculations being not controlled approximations in that error bars can 
not be put on the predictions. Although there are some attempts to resolve these problem in terms of resumming 
the momentum expansion at vector meson scale ||, the resummation was not yet used in phenomenological studies 
of vector mesons so far. The purposes of this present paper are to study chiral expansion in single-vector-meson 
processes systematically in the framework of chiral quark model, and to overcome the difficulties mentioned in above. 
Our studies will be not only resumming momentum expansion of vector mesons, but also beyond the leading order of 
A^T 1 expansion. 

The convergence of the chiral expansion is a fundamental requirement for any chiral model. In this paper, a broad 
class of chiral quark model will be checked by this requirement. The simplest version of the chiral quark model 
is in ref. ||, that was originated by Weinberg ||, and has been extended to some different versions 0-|lO|. The 
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constituent quark model in ref. [pi is different from other versions in refs. [0-|lO[. The massless pseudoscalar mesons 
are treated as fundamental dynamical field degrees of freedom in the former, but are treated as composite fields of 
quarks in the latter. In section 5 of this paper, we will show that the chiral quark models in refs. J^lOt can not 
yield convergent chiral expansion in single-vector-meson processes, and the effective lagrangian derived from these 
chiral quark models can not give unitarity S-matrix for vector meson physics. Therefore, in this paper, our studies 
will be within the framework of a chiral constituent quarks model proposed by Manohar and Georgi || (quote as MG 
model hereafter). It is the simplest and practical model for achieving our goal. In view of this model, in energy 
region between chiral symmetry spontaneously breaking(CSSB) scale(A x ~ 2^/^ ~ 1.16GeV) and the confinement 
scale(AQcr> ~ 0.1 — 0.3GeV), the dynamical field degrees of freedom are constituent quarks (quasi-particle of quarks), 
gluons and Goldstone bosons associated with CSSB(it is well-known that these Goldstone bosons correspond to lowest 
pseudoscalar octet). In this quasiparticle description, the effective coupling between gluon and quarks is small and 
the important interaction is the coupling between quarks and Goldstone bosons. The effective action describing 
interaction of Goldstone bosons can be obtained via integrating out freedoms of quark fields. 

The MG model and its extension have been studied continually during the last fifteen years and got great success 
in different aspects of phenomenological predictions in hadron physics . Low energy limit of the model has been 
examined successfully Qjl^]. In this present paper, the MG model will be extended to include lowest vector meson 
resonances (quote as extended Manohar-Georgi model or EMG model hereafter) in terms of WCCWZ realization for 
spin-1 mesons p^JT^ ]. The advantages of approach of the MG model are that, in principle, high order contribution 
of the chiral expansion can be computed completely, and only fewer free parameters are required. So far, however, 
there are still no any systematic studies beyond the lowest order of the chiral expansion due to technical difficulties. 
After freedoms of fermions are integrated out, traditionally, the determinant of Dirac operator are regularized by 
some standard methods, such as Schwinger proper time method Jl5| , or heat kernel method Jl6[ |, and are expanded 
in powers of momentum of mesons. In terms of this method, the effective action can keep chiral symmetry explicitly 
order by order. However, it is very difficult to capture arbitrary order contribution in chiral expansion by means of the 
above method. In this present paper, we will propose another equivalent method to evaluate the effective action via 
calculating one-loop diagrams of constituent quarks directly. We call it as proper vertex expansion method. By this 
method, the effective action is expanded in numbers of external vertices instead of external momentum. The effective 
action with n external vertices is just connective n-point Green function of EMG model in presence of external fields. 
It includes all order contribution of the momentum expansion. Thus in terms of this method, the chiral theory at 
vector meson energy scale can be studies up to all orders of the momentum expansion. 

Another challenging subject is to study vector meson physics beyond large N c limit. In the other words, how 
can we calculate meson loops systematically and consistently? At vector meson energy scale, it can be expected 
that the contributions from pseudoscalar meson loops play important and highly nontrivial role. This is due to 
the following reasons: 1) Empirically, for vector meson physics, the contribution from pseudoscalar meson loops is 
about T p : m p ~ 20%. The conclusion is obtained because the width of p-meson is generated dynamically by pion 
loops and T p : m p is suppressed by A^T 1 expansion according to standard N~ x argument |0. It indicates that 
the contribution from pseudoscalar meson loops is not small so that it can not be omitted simply. 2) Since both 
pseudoscalar mesons and constituent quarks are fundamental dynamical field degrees of freedoms at this energy scale, 
theoretically, loop effects of pseudoscalar mesons should be also considered as well when effective action is derived by 
loop effects of constituent quarks. 3) For physical value N c = 3, the A^T 1 expansion is not convergent very rapidly. 
Due to the above reasons, in this paper we will also study contribution from pseudoscalar meson loops systematically. 
The calculations on meson loops will suffer ultraviolet (UV) divergent difficulty. At vector meson energy scale, loop 
effects of pseudoscalar mesons cause (f> — u mixing which will destroy OZI rule if this contribution is not to be an 
appropriate value. In the other words, it can insure OZI rule for <p to choose an UV cut-off appropriately. Thus, very 
fortunately, the OZI rule for 4> provides a natural way to determine all UV-dependent contributions in the meson loop 
calculations. 

The unitarity of S-matrix is a fundamental requirement for any quantum field theory. It is well-known that the 
unitarity of S-matrix of an effective theory can not be insured in full energy region, but only be insured in energy 
region below an energy cut-off. The consistency of this effective theory requires masses of all resonances should 
be smaller than this cut-off. In the other words, it indicates that the cut-off Aemg hi EMG model should satisify 
Aemg > niuj for SU(2) and Aemg > m 4> for SU(3) at least. In general, a low energy effective field theory of QCD is a 
well-defined perturbative theory in N~ expansion. Thus the unitarity of S-matrix of low energy effective field theory 
has to be satisfied in order by order in A^" 1 expansion. For EMG model, the unitarity of S-matrix will be examined 
explicitly in this paper. We will see that the meson-loop effects play an essential role for meeting the requirement of 
the unitarity. 

The phenomenology of vector meson physics is very rich. For example, the vector meson dominance (VMD) and 
universal coupling for vector mesons, which are good approximation for vector meson physics p9| , should be satisfied 
in the leading order of the chiral expansion. In this present paper, we will achieve various phenomenological results: 
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l)The decays for p — ► tttt, p — > e + e". 2) Dynamical width for p-meson and Breit-Wigner formula for propagator of 
resonance. 3) p° — u> mixing and u> — > tt + tt~ decay. The mass difference between p and ui are also predicted. We 
will show that, for resonance particles with large width(e.g., p), the mass parameter in its propagator is not just the 
its physics mass. 4)Pion form factor and 1=1, 1 = 1 pion-pion phase shift. Some nontrivial phase in pion form 
factor are revealed too. 5) Individual light current quark masses at vector meson energy scale. The quark masses are 
model-dependent and are necessary to extend our study to K* and 4> mesons. 6)The anomalous decays p ± — ► 77r ± 
and u) — * 77r. All above results resum all information of the momentum expansion of vector mesons and are up to the 
next to leading order of iV" 1 expansion. 

The paper is organized as follows. In sect. 2 the Manohar-Georgi model and its extension are reviewed. In sect. 
3, the relevant effective action is derived at the leading order of expansion. This effective action includes all 
order contribution of the chiral expansion in powers of vector meson four-momentum square. The calculation on 
pseudoscalar meson loops will be included in the following sections. In sect. 4, decays for p — > tttt and p° —> e + e~ 
are studied. In sect. 5, the unitarity of S-matrix is examined explicitly and Breit-Wigner formula for p-propagator is 
obtained. In sect. 6, the p° — to mixing and lu —> tttt decay are studied. The isospin breaking parameter — m u is 
predicted at energy scale p, = m p . In sect. 7, the anomalous decays for p ± — > 77r ± and ui — > r yrr are predicted. In sect. 
8, pion form factor, 1=1, I = 1 tt — tt phase shift and p° — u> mass difference are studied. The light quark masses are 
predicted in sect. 9 and a brief summary is in sect. 10. 



II. MANOHAR-GEORGI MODEL AND ITS EXTENSION 



For interpreting physics below CSSB scale, Manohar and Georgi provides a QCD-inspired description on the simple 
constituent quark model. At chiral limit, it is parameterized by the following SU(3) V invariant chiral constituent 
quark lagrangian 

C x = i^{p + F + gjfon)!) - + ^Tr^V^C/V^}. (2.1) 

Here Trj- denotes trace in SU(3) flavour space, ip = (u,d,s) are constituent quark fields, g A = 0.75 can be fitted by 
beta decay of neutron. The A M and are defined as follows, 

Am - ^U f (^ - ir M )£ - ^(d, - il^}, 

IV = \{^(d„ - irjt + ^ - « f }, (2.2) 
and covariant derivative are defined as follows 

V M t7 = d^U - ir^U + iUl„ = 2£A M £, 

V M C/ f = d,U^ - a M U* + iU^r, = -2£ + A£ f , (2.3) 

where = v p + a M and = v, — a M are linear combinations of external vector field and axial-vector field 
a^, £ associates with non-linear realization of spontaneously broken global chiral symmetry G = SU(3)l x SU(3)r 
introduced by Weinberg jl3| , 

<te£(*)fc + (3) = g L ,g R eG, h{$) e H = su(3) v . (2.4) 

Explicit form of £($) is usual taken 

£($) = exp {i\ a <$> a {x)/2}, [/($) = (2.5) 

where the Goldstone boson $ a are treated as pseudoscalar meson octet. The constituent quark fields transform as 
matter fields of SU(3) V , 

— 4> — >^h\<$>). (2.6) 
is SU(3) V invariant field gradients and r p transforms as field connection of SU(3) V 

a m — ► ft($)A M ^(*), r M — > *i($)i>t(*) + h($>)d,h\<s>). (2.7) 
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Thus the lagrangian( 2.1) is invariant under G g i bai x Gi oca i- 

For achieving the purposes of this present paper, MG model must be extended to including lowest vector meson 
resonance and beyond chiral limit. The light quark mass matrix M. = diag{m u , mj, m s } is usually included into 
external spin-0 fields, i.e., x = s + Wi where s — s C xt + -M, s cx t and p are scalar and pseudoscalar external fields 
respectively. The chiral transformation for x is 



X 



9 R xg\- 



Thus x an d x^ together with £ and $ can form SU(3)y invariant scalar source £^x£t 
(£ x£ — £,X*0l5- Then current quark mass dependent term is written 



(2.8) 

£x £ pseudoscalar source 



(2.9) 



Eq. (2.9) will return to standard quark mass term of QCD lagrangian, ipMip, in absence of pseudoscalar mesons at 
high energy for arbitrary n. It means that the symmetry and some underlying constrains of QCD can not fixed the 
coupling between pseudoscalar mesons and constituent quarks. Hence k is treated as an initial parameter of the model 
and will be fitted phenomenologically. 

From the vi ewp oint of chiral symmetry only, an alternative scheme for incorporating vector mesons was suggested 
by Weinberg and developed by Callan, Coleman et. al jl4|. In this treatment, vector meson resonances 
transform homogeneously under SU(3)y, 



V M - h(<S>)V^(<t>), 



(2.10) 



where 



V ll (x) = X-V ll = V2 



( !§L 
1 V2 



V2 



\ 



K*r 



K*° 



K*+ \ 



(2.11) 
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Then extended Manohar-Georgi model(EMG model) is parametered by the following SU(3) V invariant lagrangian 



C x = iq 



(P + f + gjfo - if)q - mqq - ]-q(^X^ + - ~ ZxHH 
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(2.12) 



We can see that there are five initial parameters g A1 to, k, F and too in EMG model. These parameters can not 
be determined by symmetry and should be fitted by experiment. The effective lagrangian describing interaction of 
vector meson resonances will be generated via loop effects of constituent quarks. 

It should carefully distinguish EMG model from ENJL-like models |||],|2(J which seem to be similar to EMG model. 
The difference is that there are kinetic term of pseudoscalar mesons in EMG model, but this term is absent in ENJL- 
like models. It makes there is a extra constrain on constituent quark mass m in ENJL-like models, and this constrain 
requires constituent quark mass m around 300MeV. However, in sect. V we will point out that this value of m can 
not insure convergence of chiral expansion in single-vector-meson processes. 



III. EFFECTIVE ACTION AND PROPER VERTEX EXPANSION 



A. Proper vertex expansion 



The essential subject of nonperturbative QCD is to obtain a low energy effective action of light hadrons, which 
belongs to the leading order of A*"" 1 expansion. In chiral quark model, this effective action is generated via loop 
effects of quarks, i.e., integrating over degrees of freedom of quark. This path integral can be performed explicitly. 
The determiant of fermion obtained by path integral are usually regularized by heat kernel method. In this method, the 
effective action is expanded in powers of momentum of mesons. It is very difficult to use this method to calculate high 
order contribution of momentum expansion. In this present paper, we will derive the effective action via computing 
the feymann diagrams of quark loops directly instead of performing path integral. 



4 



We start with constituent quark lagrangian (2.12), and define vector auxiliary field VHa = 0, 1, • • ■ , 8), axial- vector 
auxiliary field A" scalar auxiliary field S a and pseudoscalar auxiliary field P a as follows 



S a = ^Tr f {\ a (£}x£} + &Qh P a = jTr f {X a ^ X ? ~ &Z)h (3-1) 

where A 1 , • • • , A 8 are SU(3) Gell-Mann matrices and A = y|. Then in lagrangian ( [2.12 ), the terms associating with 
constituent quark fields can be rewritten as follow 

£* = q(ip - m)q + V?q\ a >fq + ig A A-qX a ^ l5 q - S a qX a q - P a qX a l5 q. (3.2) 

The effective action describing meson interaction can be obtained via loop effects of constituent quarks 

e iS -« = <0\T q e l f d4xC * {x) \0> 

= Yx i f rf4pi ^^ " "^? ] ^ n ^ 1 '" ■' p ™^ 4 ^ 1 ~ P2 Vn "> 

= ffl l( 0) + f> / ^-^^{pu-^n-!), (3.3) 



n=2 



where T q is time-order product of constituent quark fields, C l x is interaction part of lagrangian ( |3.2j ), n n (pi, ■ ■ ■ ,p n ) 
is one- loop effects of constituent quarks with n auxiliary fields, pi,P2, •• • ,Pn are four-momentas of n auxiliary fields 
respectively and 



n„(pi, • • • ,p n -i) = / p„n n (pi, • • • ,p„)<5 (pi - p 2 Pn)- (3-4) 

To get rid of all disconnected diagrams, we have 



n=l 

5i=n!(0), (3.5) 
f d 4 Pl d 4 p„_i 

Obviously, in eq. ( |3.5| ) the effective action S e g is expanded in powers of number of external vertex and is expressed as 
integral over external momentum. Hereafter we will call this method as proper vertex method, and call S n as rt-point 
effective action. 



B. Relevant effective action 



In terms of proper vertex method, the effective action can be obtained, and in principle it can include all information 
on resumming momentum expansion. Practically, however, it is tedious to express the effective action completely. In 
this section, we only derive those effective action relating to phenomenological studies on p and ui physics in next 
sections. For simplifying the expression of the effective action, some good approximation will be used. Since we work 
in p-meson energy scale, the expansion of light quark masses still works very well. We can set m u = = in 
p physics and take leading order of — m u expansion in isospin breaking lo physics. This also implies that the 
soft-pion theorem for pseudoscalar mesons is available here. 

The one-point effective action S\ is just contribution of tadpole diagram of constituent quark, 

Si = J d 4 x^-B Tr f { X U^ + X f U}, (3.6) 
where Bq is a constant which absorb the quadratic divergence from quark loop integral 
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(4tt) b / 2 



(3.7) 



After renormalizing the kinetic term of pseudoscalar mesons, the two-point effective action reads 
S 2 = ^|y d^xTrfiV^UV^} + ^rnl J d 4 xTr f {V p (x)V^x)} 

" ' <J " 2) -(W 2 -P*P»)Trf{(V»(p) + iT»(p)){V»{-p) + iT»(-p))} + 



(27r) 4 4 w ' x « / *' riurvj- 1 - • j 1 \' \f) ' \t>jj\" \ f) > " v -f^yyj 1 ■••) (3-8) 
where denotes those terms are independent of the following phenomcnological studies of this present paper, and 

«i(p 2 ) = 5 2 - -| / -x 1 - a:) In 1 - -i 2-^- , 

0(p) = / d 4 xe ip - x O(x), O = V, T, A, • • • . (3.9) 



The constant g in eq. ( |3.9| ) is an universal coupling constant, which absorbs the logarithmic divergence from quark 
loop integral, 



Nr. 



_„- ^£ fifl^/2 rr2 _ ^\ 



(3.10) 



It is obvious that the second term of S2 includes all order information of momentum expansion. In coordinate 
space, it can be expanded in powers of number of derivative, 



-Tr f {% v Vn 



Nr. 



240tt 2 



2 Tr f {d p V^d<>V^} 



N r 



2240tt 2 



m- 4 Tr f {d p d„V llv d p & T V'"'} 



with 



%u = d*% - d v Va 



V ll = V ll +iT ll . 



(3.11) 



(3.12) 



In next section we will show that the high order derivative terms yield very important contribution at vector meson 
energy scale, i.e., for p 2 ~ m 2 . 

The three-point effective action includes two parts. One describes reactions with normal parity and another one 
describes reactions with abnormal parity. For the purpose of this present paper, the three-point effective action with 
normal parity can be obtained via standard calculation of fcynman diagram(fig.l), 





S. 



(normal) 



FIG. 1. The triangle diagrams of constituent quarks. 



(^ 2 - 



N c 
16tt 2 



)mg\ d^xTrfiV^UV^ixrf + Utf)} 



+ 



d 4 x 

Nr 



d 4 P 



ip x {--^-a 2 (p 2 )p" Tr, [V v (p) [A„ (x) , A„ (x)}} 



(2tt) 4 l 2 

a 3 (p 2 ) (g^ P 2 - P»Pv)Tr f [V" (p) V" (x) € + th^ )] 



247r 2 m 

^*9 A «4 {p 2 )Tr f [W (p) , A M (x)] - ftft)] 
AT, 



167r 2 m 



rzg A a 5 (p*)( gfiV p 2 -p IM p u )Tr f [V^p),A v (x)}(^x^ -£x 1 0}} + 



(3.13) 



G 



where the form factor cti(i — 2, 3, 4, 5) are defined as follows, 



N f 1 f 1 m 2 f(v 2 ) 

a 2 ( P 2 ) =9>-^j o d X I dy • ,(1 - X y)[l + + ln(l - ™2)], 



Oi3(p 2 ) = dx 
Jo 



2 



6x(l — x)m 
m 2 — x(l — x)p 2 ' 

1 <" ! f(p 2 ), , xy(l-x)p 2 



with 



a 4 (p 2 ) = 6*V - ^ dx dy{2xHl - ^_Z) + <;-,. N ) 

, 2N f 1 , f 1 , 2x{l-x)m 2 
a 5 (p)=] Q dxJ Q dy m2 _ f(p2) , 



f(p 2 )=x(l-x)(l-y)p 2 . (3.15) 



Moreover, the following three-point effective action with abnormal parity is related to the purpose of this present 
paper, 

^ = ^9^ J ( d 2 4 J ) > 6 (gLg|)gi^ / {K(gi)^(g 2 )AMgi +92)} + -, (3.16) 

where 

«e(g?,g£) = f dx ■ x [ dy{\ - (x + y)(l - x)(l - y)^±A}~^ (3.17) 
Jo Jo zm 

It has pointed out that, the abnormal part of these quark models including spin-1 meson usually do not pose the 
same symmetry as the QCD anomalous Ward identities |2l| , and the authors of ref. [^2| have provided a method 
to subtract the extra terms. However, in this present paper we only focus our attention on anomalous p — ► jtt and 
lu — > 77r decays. Thus we have notation A^fc) — * kpir(k) in eq. (3.16). Using this notation we have 

e^ al3 Tr f {q ltl V„( qi )r a (q 2 )A fj (k) + q^V^T a {q 2 )Ap{k)} 
- e^eTr } {-q^V v {qi)Ql a {q2)^{k) + qi M Qrf v (qi)V a (q2)A (k)} (3.18) 

be invariant under infinite transformation 

Ol^O% lfi ^^ + tdfl e°, (3.19) 

where e c , (c = 1, 2, ...8) and e° are infinite parameters. Therefore, here no extra terms need to be subtracted. 

The four-point effective action also includes both of normal parity part and abnormal parity part. Due to soft-pion 
theorem, the normal parity part relating to this paper reads 



Nc 2 f , 4 f A 



X \ (2^ e ~ iq ' X{9 ^ ~ 9 »° qv) {^{q 2 )Tr } {{V^q), (£^ + ^jA^A^)} 



8n 2 m aA J J (2tt) 4 

-a s (q 2 )Tr f {V»(q)A»(xmH + tW^i*)}) + (3-20) 



where 



a r(Q 2 ) = / dxi ■ x\ j dx 2 (l — x 2 ) 
Jo Jo 

ots{q 2 ) = / dx 1 -x\ I dx 2 (l-x 2 ) 



3 - 2x\x 2 (l + 2xi)(l - x 2 )q 2 /m 2 

[I - x\x 2 (\ - x 2 )q 2 / m 2 ] 2 
4(1 - Xl )[3 - 4x 2 x 2 (l - x 2 )q 2 /m 2 } 
[1 — x 2 x 2 (l — x 2 )q 2 1 m 2 ] 2 
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For abnormal parity part, we only focus on photon to three-pseudoscalar interaction in this paper. Thus using 
soft-pion theorem, we have 

^ r3$ = / S " ' ^i^ VaP ^^Ml2 + 93 + q,)Tr f {Q<S>( q2 )<S>( q3 )<P( qi )}, (3.21) 

where $ is pseudoscalar meson fields, is photon field, Q = diag{2/3, — 1/3, — 1/3} is charge operator of quark 
fields. 

So far, all meson fields in above effective actions are still non-physical. The physical pseudoscalar meson fields 
can be obtained via field rescaling P — > 2F~ X P (here P denotes pseudoscalar mesons, and we ignore the difference 
between F v and Fk since we focus on p and to physics only). For defining physical vector meson fields, let us pay 
attention to kinetic term of vector mesons, which can be obtained from eq. (3. 



412 = -~S 2 2>/{(^K - d v Vj{d»V - d»V»)} + \g 2 m 2 v Tr f {V^}. (3.22) 

Here 

o ml N r n f 1 x(l — x)m 2 

m 2 =^ + ^™J <fe^(l-x)ln(l-- J—*-) (3.23) 

g z 7i"V v Jo m 

is a "common" mass parameter of vector meso ns. It implies that the masses of light vector meson octet are degenerated 
at chiral limit and large N c limit. From eq. ( |3.22 ) we can see that the physical vector meson fields can be obtained 
via field rescaling — > g^ 1 V fl . 



C. Vector meson dominant and KSRF sum rules 



The direct coupling between photon and vector meson resonances is also yielded by the effects of quark loops. 
Therefore, when vector meson resonances are treated as bound states of constituent quarks, vector meson dominant 
will be yielded naturally instead of input. At the leading order of large N c expansion the VMD vertex reads from 
eq. (|J) 

£vmd = -\e J ^e-^f$(q 2 )(g^ q 2 - q , q „)A»(x)Tr f {QV»(q)}. (3.24) 

where fpj\ q 2 ) = g~ 1 ai( q 2 ). At the leading order of momentum expansion, i.e., fjry(q 2 — 0), the above equation is 
just the expression of VMD proposed by Sakurai In addition, at the leading order of large N c ex pansio n the 

V — vertex (where V denote vector mesons and $ denotes pseudoscalar mesons) reads from eqs. (HJ) and ( §lf ) 

£v** =~J ^e-^f$l( q 2 )(g^q 2 q^)Tr{V^q)mx), V*(x)]}, (3.25) 

where 

f%W) = -^lai(l 2 )-9> 2 ( q 2 )} (3.26) 

It is well known that the KSRF (I) sum rule §|| 

f P ,( q 2 ) = f P Mq 2 )K (3-27) 
is expected to be available at the leading order of momentum expansion. At this order we have 

/(°V = 0) = . 9 , 



»y 2 = ')) = ^{s a + (^-/r)}- I*-*) 



We can see that KSRF(I) sum rule is satisfied exactly for g — y Therefore, 9 — y (especially, g = ir 1 for 

N c = 3) is favorite choice for the universal constant of the model. 
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D. Low energy limit 



It is well known that, at very low energy, ChPT is a rigorous consequence of the symmetry pattern of QCD and 
its spontaneous breaking. So that the low energy limit of those model including meson resonances must match with 
ChPT. The low energy limit of EMG model can be obtained via integrating over vector meson resonances. It means 
that, at very low energy, the dynamics of vector mesons are replaced by pseudoscalar meson fields. Since there are no 
interaction of vector mesons at 0(p 2 ), at very low energy, the equation of motion SC/SV^ = yields classics solution 
for vector mesons are follow 

Vu = —7T x 0(p 3 )terms, (3.29) 

V 

where p is momentum of pseudoscalar at very low energy. Therefore, in effective action S n , the terms involving 
vector meson resonances are 0(p 6 ) at very low energy and do not contribute to 0(p 4 ) low energy coupling constants, 
Li(i = 1, 2, 10). The low energy coupling constants Li(i ^ 7, 8) yielded by EMG model can be directly obtained 
from effective actions in subsection 3.2, 

T _ 1 T N c _ N c N c 4 

Ll _ 2^ ~ 384^*' Ls ~ "64^ + 192^ ffA ' 

La = L 6 = 0, L 5 = T^grffi (3.30) 

?2 ™ „2 



N n r Nr. N c 



— ttt^j L w — — u + " 2 3a- 

In fact, the above expression on Li have been obtained in some previous refs. (besides of L$). 

The constants L7 has been known to get dominant contribution from 770 H and this contribution is suppressed by 
l/N c . If we ignore the 77 — 77' mixing, we have 



f 2 



128m2, ' 



(3.31) 



Thus the six free parameters, g(it has been fitted by KSRF sum rule), g A (it has been fitted by n — > pe v e decay), 

Bq, k, m and determine all ten low energy coupling constants of ChPT. It reflects the dynamics constrains between 

2 

those low energy coupling constants. Moreover, if we take m u + rrid — lOMeV, we can obtain Bq = — m + m — 2GeV. 
Then experimental values of Lg constrains constituent quark mass m ~ 480MeV (a more detail fit is in sect. IX). 
The numerical results for these low energy constants are in table 1. 





Li 


L 2 


L 3 


U 


L 5 


U 


L 7 


L 8 


L 9 


L\o 


ChPT 


0.52 ±0.23 


0.72 ± 0.24 


-2.7 ± 1.0 


-0.3 ±0.5 


1.4 ±0.5 


-0.2 ±0.3 


-0.26 ±0.15 


0.47 ±0.18 


6.9 ±0.7 


-5.2 ±0.3 


EMG 


0.79 


1.58 


-4.25 





1.26 





(-0.4±0.1) f,) 


0.47 a) 


6.33 


-4.55 



TABLE I. Li in units of 10 3 , fj, = m p . The experimental data is from ref.[24]. a)input. b)contribution 
from gluon anomaly. 
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E. Chiral expansion at m p -scale 



Let us calculate the on-shell decay p — > irir and p° — ► e+e~ under the leading order of momentum expansion 
and under including all order information of vector meson momentum expansion respectively. Numerically, using 
g A = 0.75, g = ir^ 1 and m = 480MeV, we have 

/#V = 0) = ^ = °- 318 > 4°l(9 2 = 0) = (0-93 x 10- 5 )MeV- 2 , 

f$(q 2 = m 2 p ) = 0.439, f^(q 2 = m 2 p ) = (1.14 x 10- 5 )MeV" 2 . (3.32) 

Then the decay widths are T(p — > 7T7t) = 125MeV and — » e + e~) = 4.35MeV at leading order of momentum 
expansion, and T(p — > 7T7t) = 188MeV and r(p° — > e+e~) = 8.29MeV when high order contributions of momentum 
expansion are included. It implies that the high order contributions of momentum expansion are very important at 
vector meson energy scale. Furthermore, it indicates that the chiral expansion converge slowly at vector meson energy 
scale. This is an important feature of chiral theory including vector meson resonances. 

However, neither the results at the leading order of momentum expansion nor the results including high order 
contributions of vector meson momentum expansion match with experimental data, r cxp (p — > itit) = 150McV and 
Tcxplp — > e + e~) = (6.77 ± 0.32)MeV. The reason is that, so far, the contributions from pseudoscalar meson loops 
are still not included. We have pointed out in introduction that these -/V" 1 suppression effects yield about 20% 
contribution. Therefore, in the following sections, we will calculate the effects of pseudoscalar meson loops for various 
physical processes. 

IV. ONE-LOOP CORRECTION TO p DECAYS 

In this section we will calculate pseudoscalar one-loop correction to p — > nir and p° — > e+e~ decays and provide a 
complete prediction on these reactions. Since m 2 <C m 2 K < m 2 , we will treat pion as massless particle but m K ^ in 
our calculation. In addition, we will omit the difference between m n8 and m K , and between F„ and F K , since these 
differences are doubly suppressed by light quark mass expansion and iV" 1 expansion. 

A. Correction to p — > nn decay 

To the next to leading order of iV c _1 expansion, there are three kinds of loop diagrams of pseudoscalar mesons need 
to be calculated (fig. 2). Recalling m 2 = but m 2 g = m 2 K ^ 0, we can see that only K and i] 8 mesons can yield 
non-zero contribution of tadpole diagram in fig. 2-a), since in dimensional regularization, J d D kk 2n = 0(n > — 1). 
Due to the same formula, we have 

f k 2s 

/ d D k— cx m 2s x quadratic divergence, (s > 1). 

J k 2 +m^ + ie v 

Thus in this present paper we will ignore all contribution from quartic divergence or higher order divergences. Fur- 
thermore, in the following calculation we can find that the lowest order divergence is quadratic divergence. In means 
that the on-shell condition for pseudoscalar mesons and soft-pion theorem in tree level are still available even in 
calculation on pseudoscalar meson loops. In addition, when we calculate two-point diagram of pion, we must include 
all chain- like diagrams in fig. 2-c). Meanwhile, when we calculate two-point diagram of K or rjg,, we only need to 
calculate one-loop diagram in fig. 2-b). The reason is that the two-point diagram of pion generates a large imaginary 
part of T-matrix for p and lo physics, but K or r/s do not. If we only calculate one-loop diagram in fig. 2-b), a large 
error bar for phase factor will be yielded and this error will affect the study on pion-pion phase shift in section 7. 
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K,T] S K 
P C ) n _Q_ 





(a) (b) 



+ K + A + 




(c) 

FIG. 2. One-loop diagrams for p — » 7T7t decay. a)Tadpole diagram of if or 778. 
b)Two-point diagram of K. c)Chain-like approximation of pion. 

It is convenient to evaluate pseudoscalar meson loops in terms of background field method. To expand pseudoscalar 
meson field around their classic solution 

U(x) = £»e^», Um = ?(*), (4.1) 

where background field f ($) is solution of classic motion of pseudoscalar mesons, 6£/5U(x) — 0, ip(x) is quantum 
fluctuation fields around this classic solution. 



Let us first calculate the tadpole diagram (fig. 2-a). Inserting eq. (4.1) into effective action S2 and S3 and to retain 



terms including two $ and two quantum fields <p, the p — 4$ vertex reads, 

£ p~4* = -^5 / ^e-^f^{ q 2 )Tr f {p^M{d^d^^ + V d^d v $<p - d^d^tp, ^})} + • • • . (4.2) 

where q are four- momentum of p, p^ v {q) = q^Puio) ~ ^vPni*}), and " • • •" denotes those p— 4$ terms that only generate 
zero contribution, e.g., Tr f{p fiv (q)d^(p^d 1 ' <p&} . 

In terms of the completeness relation of generators X a (a = 1, 2, • • • , N 2 — 1) of SU(N) group 

2 

Tr{\ a A\ a B) = — -Tr(AB) + 2TrATrB, 

Tr{X a A)Tr{X a B) = 2Tr(AB) - ^TrATrB, (4.3) 

we have 

Tr f {p^(q){d^^X a d u ^X a + X a d t *$>d v <$>X a - d^d v ${\ a , X a })} = -4NTr f {pn V (q)d»<f>d v <f>}. (4.4) 

From viewpoint of group theory, the contribution from K and r/g, is just one from pseudoscalar mesons in SU (3) / SU (2) 
sector. Then tadpole diagram correction to p — > is 



= / ^e- i «-/W 7r ( g 2 )(^ g 2 -g /1 g„)Tr / {^(9)[$W,^$(a ; )]}, (4.5) 

where A is a constant absorbing quadratic divergence from meson loops 

AH^V 2 F(l-f). (4.6) 
m z 2 

In the following we will calculate the two-poi nt di agram(TPD) of pseudoscalar mesons. The leading order 
vertex of N^ 1 expansion has been given in eq. ( [3.25| ). Here we need to further obtain explicit 4-$ vertex from the 
effective actions of section 3. Inserting eq. (4.1) into S2, s| nor:ma1 ' and g l ^ lolmal '> anc j retaining terms with two <I> and 
two ip fields, we have 
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(2ny 

+ •••, (4.7) 

with 

1 



^ = 77^2^^]- (4-8) 



Here " • • •" denotes that those 4-$ terms only yield zero contribution due to transveral operator g^q 2 — q^qv in 
V — vertex, such as Tr/{[9 M <I>, <p] [9 M< &, ip]}. Then two-point diagram correction to p<I>$ vertex is 

= / /™ rf «( 2 ') 4 /i;l(? 2 )fc,9 2 - 9„fc)fo*> + jfj(s 2 - " 4«9«)]rr/{p"(4)[A", A»]) 

X r,- J{ [. M , y . w „A- , A1> / g - 5 - - . (4.9) 



Using completeness relation of generators of SU(N) group (4.3), we have 

7>{A[A a , A b ]}Tr{ J B[A a , A b ]} = -8NTr(AB) + 8TrATrB. (4.10) 
Then taking N = 2, we obtain SU(2) (pion meson) correction as follow 



where 



with 



£$»g[Stf(2)] = \J ^e-^'fWtf^tf)^*? - q,q„)Tr{V»(q){<S>(x),d^(x)}}, (4.11) 

+ -Ur. 9 (-l)c% 2 - 4m*)(l - -^#) 3 / 2 }, (4.12) 
6 g z 

Arg(-1) = (l + 2fe)7r, fe = 0,±1,±2,- • •, 



where 



In eq. (4.12), the phase factor (1 — 4m 2 /g 2 ) 3 / 2 denotes th e pio n mass correction, which is yielded by unitarity of 
S'-matrix(see section 5). In addition, taking N — 3 — 2 in eq. ( 4.10 ), we can obtain SU(3)/SU(2) (if -meson) correction 

C$Z[SU(3)/SU{2)] = --J ^e-^f^{q^ K {q 2 ){g, v q 2 q,q u )Tr{V^q)Mx),d^(x)}}, (4.14) 

„ , on m i r 3o 2 . n »1 ,,„, o 2 , f 1 , . x(l — x)q 2 , , , x(l ~ x)q 2 ^ 

= g^[l + ill - ^)]{A(1 - ^) + / ^(1 - -4^-) Ml - ^F^)i- ( 4 - 15 ) 

Comparing eq. (4.11) with eq. ( 3.25| ), we can see that every two-point loop in chain-like approximation (fig. 2-c) 
contributes a factor — E^g 2 ). Therefore, to sum over all diagrams in fig. 2-c), we obtain 

££T = —J J^ e ~ iq ' X YT^m [9 ^ q2 ~ W*>)Tr{V»(q)Mx),tr*(x)]}. (4.16) 



Eq. (4.16) together with eqs. ( |4.5| ) and (4.13) lead to "complete" p — tttt vertex to the next to leading order of iV c 1 
expansion, 

* f dA<1 r- 1 
(27T) 4 



41 = -I / Ve-^ffrKq 2 )^ 2 - q,q,)Tr{p^q)[7:(x),d^(x)]}, (4.17) 



where 
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B. Correction to p — > e + e decay 

To the next to leading order of N^ 1 expansion, we need consider three kinds of loop diagrams of pseudoscalar 
mesons correction to VMD vertex (fig. 3). 





b) 



VW + — I^Wv + 

c) 

FIG. 3. One-loop diagrams for p° — 7 mixing. a)Tadpole diagram of K or 77s . 
b)Two-point diagram of K. c)Chain-like approximation of pion. 



We first evaluate the tadpole contribution in fig.3-a). Inserting eq. ( }4.1| ) into S2, we obtain p^/ip^p vertex as follow 
£ P7W = 4 e / 4k e ' iq ' X f^^)i9^q 2 - q»q v W{x)Tr s {p»{q){ l -Q{ V ^} - <pQ<p)}. (4.19) 



F 2 J (27r) 4 J m^ >\»^ v-v-'jir wv 2 

Then the tadpole correction of K or 7y 8 mesons to VMD vertex is 

4Ld = l^k e J ^^f^il^q 2 qtt q„)A»(x)Tr f {Qp»(q)}. (4.20) 
For calculating contribution of fig. 3-b) and fig. 3-c), the 7 — ipip vertex is necessary 

C-rw = ~e J (0e"^[l + f/jSlfa 2 )^ - ^M^Tr^eM*). d^(x)]}. (4.21) 
Calculations on fig. 3-b) and fig. 3-c) are similar to one on fig. 2. The final result on "complete" VMD vertex is 
47 = -\\ ^ e ~^/£V) W ~ q»qu)A»(x)Tr f {Qf{q)}, (4.22) 



where 



with 



/&V) = 4°V)(i - - f£l(« %lf2 2 ) " ^ l[q2% (4 ' 23) 



x(l - x)q 2 



^{q 2 ) = ^[1 + f <? 2 /<°l(<7 2 )]{~ + ^ *fa • *(1 ~ *) In 
+ ^r ff (-l)% 2 -4m^)(l-M)3/2 }j 

^) = g[l + |^)HA(1 - ^) + jf -(1 - ^) Ml - ^) } . (4.24) 
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C. Cancellation of divergence 



From the above calculations we can find that only quadratic divergence appears in one-loop contribution of pseu- 
doscalar mesons. Since the present model is a non-renormalizable effective theory, the divergences have to be factorized, 
i.e., the parameter A has to be determined phenomenologically. 

The on-shell decay <j> — > rnr is forbidden by G parity conservation and Zweig rule. Experiment also show that 
branching ratios of this decay is very small, B{4> — > nn) = (8 _ A ) x 10~ 5 . Theoretically, this decay can occur through 
photon-exchange or iGloop(fig. 4). The latter two diagrams yield non-zero imaginary part of decay amplitude. Thus 
the real part yielded by the latter two diagrams must be very small. We can determine A due to this requirement. 
From the calculation in the above two subsection, we see that result yielded by the latter two diagrams is proportional 
to a factor 

K^-O- f dx ■ K ~ < l ~ ^)P 2 ] ln C 1 - )• ( 4 - 25 ) 

D JO m K 

Then Zweig rule requires 

{K^~rn 2 K )- R zl dx ■ K - s(l - X )P 2 ] m (1 - ~2 )T>2 )}\p-= m \ ~ °- ( 4 - 26 ) 
2 Jo m K 

Form the above equation, we obtain 

A ~ 0.54. (4.27) 






b) 



c) 



FIG. 4. Some diagrams for <j> — > mr decay. The one-loop in figure b) and c) is K-\oop. 



D. Numerical results 



The momentum-dependent widths for p — > ttt: and p — > e + e can be obtained from the "complete" vertices (4.17) 
and C22|), 



48^ 



(4.28) 



Inputting m = 480MeV, g A = 0.75, g = ir , A = 0.54 and y/q 2 = m p = 770MeV, we obtain the on-shell decay widths 



T(p -> tttt) = 146.4MeV, 
T(p^e + e-) = 6.98KeV. 



(4.29) 



These results agree with experimental data, 150MeV and (6.77 ± 0.32)KeV, excellently. It also implies that the 
one-loop diagrams of pseudoscalar mesons yield significant contribution at vector meson energy scale. 
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E. Loop correction to 7 — » tttv vertex 



Here we will further calculate loop correction of pseudoscalar mesons to 7 — * rnr vertex, which is needed by studies 
in the following section. 

The calculation of loop correcti on to 7 — > tttt vertex is similar to one to p — > rnr vertex (fig. 5). The 7$$ vertex 
and 4-$ vertex has given in eqs. (4.21) and (4.7) respectively. Here 7 to 4-$ vertex is further needed for evaluating 
the tadpole contribution in fig. 2-a), 

xTr f {Q{<$>(pd"<i>(p + Lp$d u <S>ip + <p$ipd v <& - ^$<9"${<p, <p}}. (4.30) 
Then due to conservation of electromagnetic current, the "complete" 7 — ► tttt is 

= \l ^ e ~ iq ' Xe F^)MQ)Tr f {Q[n, d^]}, (4.31) 

where F 11T7T {q 2 ) is nonresonant background part of pion form factor, 



2 



Since S^g 2 ) — Am 2 /(8tt 2 F 2 ) is 0(q 2 ), F J7T7T (q 2 — 0) = 1 is still satisfied. It means that this loop correction does not 
break U(l) e .m. gauge invariance. 



V. UNITARITY AND LARGE N c EXPANSION 



The unitarity of S-matrix, or optical theorem 

all 7 

has to be satisfied for any well-defined quantum field theory, where the 7g, Q transition amplitude, and 7 denotes 
all possible intermediate states on mass shells. It is well-known that a low energy effective meson theory should be 
a well-defined perturbative theory in -/V" 1 expansion jljj ]. Therefore, we can expand T-matrix in powers of iV^ 1 
expansion, 

00 

T = J2 Tn > T » ~ 0((^Mr n )- (5-2) 

n=o 

Then the unitarity of S'-matrix for a low energy effective meson theory can be proved order by order in powers of 
iV c _1 expansion, 

Im(7a iQ )„ = - X! {^, a ) m {T*] ) n - m . (5.3) 

all 7 m<n 

In this present paper, we will discuss the unitarity of S'-matrix of several chiral quark model at vector meson energy 
scale. In particular, as an example, we will directly examine unitarity in EMG model via forward scattering of 
p-meson. 

From effective action S2, we can see that F„ and g are both 0(\fW c ). So that every meson external line is 0(l/\fW c ) 
and every meson internal line is 0(1/N C ). In addition, from effective action S n we also see that every vertex is 0{N C ). 
Therefore, any transition amplitudes with n v vertex, n e external meson lines, rii internal meson lines and n; loops of 
mesons are order 
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Ar l-n ! -»e/2 j ( 54 ) 

where relation ni = — n v + 1 has been used. 

Using the above power counting rule, we can obtain a well-known fact that unitarity of S-matrix requires that 
on-shell transition amplitude from a vector meson to any meson states must be real at the leading order of N^ 1 
expansion. 

From eq. ( |3.26 ), we can see that, at the leading order of -/V c _1 expansion, the on-shell p — > tttt transition amplitude 



is proportional to /pji(mj). Thus the above theorem requires ai(mV) and ^(m^) defined in eqs. (3i)) and ( 3.14 ) 
are real. In the other words, it means that, in chiral quark model including vector meson resonances, the 
unitarity of S'-matrix requires constituent quark mass m > m p /2. In fact, this point is insured in EMG model 
by to ~ 480MeV which is fitted by low energy limit of the model. However, in many other version of chiral quark 
model 0HJ2^|, the constituent quark mass are constrained around 300MeV. Obviously, this value does not yield 
unitarity S'-matrix for single- vector meson processes, since all form factor defined in section 3, oti(q 2 ), are not real for 
q 2 ~ m 2 It also means that, in those model, the chir al e xpansion at vector meson energy scale is not convergent. 



In the following we will explicitly examine the eq. (|5.l|) in forward scattering of p— meson up to three-loop level of 



mesons. The examination on other processes can be p erfo rmed similarly. For a — [3 — p in eq. (5.1), < 7I =< 7T7t| is 
dominant. The for forward scattering of p-meson, eq. (^]) becomes 



2 r , m r„ = to, ™) - - <=£)*. (5.5) 



(2*)* '•' 

(c) . . _1 

The both of T PtP and can be expanded in scries of N~ expansion, 



? P , P = E( t p.p)2« 41 = E(4-Wi- ( 5 - 6 ) 



n=0 n=0 
K,% K 

p O p -p- 




(a) (b) 




O 

FIG. 5. One-loop diagrams correcting to p propagator. a)Tadpole diagram of K ot rjs. 
b)Two-point diagram of K. c)Chain-like approximation of pion. 

At the leading order, Im(7^ iP )o = is obviously satisfied. For obtaining the imaginary part of (7^ jP )2, we need to 
calculate meson loop correction in fig. 5. The calculations are similar to one in section (IV. B), and the result is 



r}-loop _ 1 f A -jg-zr >(Q) /^^r^V) 2 2 v(0) / 2 



where 



4°V) = ^~A+fdx- x(l - x) In " (1 2 X)q2 + Urg(-l)e(q 2 4m*)(l - ^!) 3 / 2 }, 
87r z 6 J m£ 6 <r 

^V) = ^{A(l-^) + r^-^^)Ml-^^)}. (5.8) 



Then taking Arg(— l) = — 7r, we obtain 



1G 



(27T)' 



-Im(T PjP ) 2 



48tt 



4to ^3/2 



(5.9) 



Noting that (/pxlr)i = /p^lr, we can see the eq. ([5~5|) is satisfied at the next to leading order. 



_D ^_ 
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a) 



b) 



c) 



FIG. 6. Two-loop and three-loop diagrams correcting to p propagator. Here loops with 
solid line denote K or -qg loops, and loops with dash line denote pion loops. It should be 
pointed out that in fig. a), b) and c), the positions of pion loop and K loop can exchange. 



If the two-loop diagrams and three-loop diagrams in fig. 6 are further considered, we can prove the eq. (5.5) is 
satisfied up to Q(N ~ 3 ). The details of calculation are omitted here. 

Finally, eq. (5.5) also tell us that the "complete" propagator of p- meson is 



AW) = 



nip + i 



q 2 T p (q 2 



(<x -^rr) term), 



(5.10) 



where T p (q 2 ) = T p ^ 7r7r (q 2 ) given in eq. ( 4.28| ). Sinc e vector meson is treated at tree level, the (ex qfiqv/m 2 ) term in 
propagator is unimportant. Then propagator (5.10) is just well-known Brcit- Wign cr formula for spin-1 resonances. 

Because the width in p-resonance (possessing a complex pole) propagator ( 5.1 0| ) is momentum-dependent, it must 
be addressed that the mass parameter rh p is not the physical mass m p = 770MeV. Let us interpret this point 
briefly. Empirically, the physical mass of resonance is defined as position of pole(real value) in relevant scattering cross 
section, or theoretically, it should be defined as real part of complex pole possessed by resonance. It is well-known 
that the width of p- resonance is generated by pion loops. For a simple VMD model, the leading order of p — ttit 



coupling is independent of q 2 



Thus one has T p VMD \q 2 ) cx 



and due to equation 



rh 2 



■ ^ p 2 

i — — q 

m„ 



0. 



(5.11) 



m 



m 2 p {l -ie + 0{e 2 )) with e = T p /m p ~ 0.19. The result yields 
784MeV. In particular, in the WCCWZ realization for vector meson fields used by this paper, 



we obtain the complex pole of p-resonance is q 
m 



p — irir coupling is proportional to q 2 at least. Hence one has T p (q 2 ) cx q 

r 



m. 



0. 



m: 



- ' q 2 at least, and complex pole equation 

(5.12) 



It yields rh„ 



+ 3e 2 = 810MeV, which poses a significant correction. 



The above discussions imply that: 1) For resonance with large width, the mass parameter in its propagator is 
different from its physical mass. The correction is proportional to the ratio of resonant width to physical mass. 2) 
The mass in the original effective lagrangian only emerges as a parameter instead of a physical quantity measured 
directly in experiment. 3) The choice of mass parameter is relied on the choice of model. But the physical quantity 
must be independent of this choice. 

Since in our result all hadronic couplings include all order information of the chiral perturbative expansion and 
one-loop effects of pseudoscalar mesons, the momentum-dependence of T p (q 2 ) is very complicate. It is difficult to 
determine rh p via the above method. The prediction on rh p will be found in section VII. 



VI. 



DECAY 



The investigation of p° — uj mixing, which is considered as the important source of charge symmetry breaking p5| in 
nuclear physics, has been an active subject pq- [R6|. The mixing amplitude for on-mass-shell of vector mesons has been 



17 



observed directly in the measurement of the pion form-factor in the time-like region from the process e + e~ — > it + it~ 
p9| , j37f . For roughly twenty years, p° — uj mixing amplitude was assumed constant or momentum independent even 
if p and u> have the space-like momenta, far from the on-shell point. Several years ago, this assumption was firstly 
questioned by Goldman et. al. |B2], and the mixing amplitude was found to be significantly momentum dependent 
within a simple quark loop model. Subsequently, various authors have argued such momentum dependence of the 
p° — uj mixing amplitude by using various approaches J33|,[54]]. In particular, the authors of ref. f34j| has pointed out 
that p° — uj mixing amplitude must vanish at q 2 — (where q 2 denotes the four-momentum square of the vector 
mesons) within a broad class of model. This point will be also examined in EMG model. 

In the most recent references ]26|-p8|], the uj — > it + tt~ decay was treated as being dominant via p-resonance exchange, 
and the direct ujtt + it~ coupling is neglected. It has been pointed out in ref. [0 that the neglect of uj "direct" coupling 
to 7r + 7r~ is not valid. It can be naturally understood since tttt can make up of vector-isovector system, whose quantum 
numbers are same to p meson. Thus in an effective lagrangian based on chiral symmetry, every p field can be replaced 
by 7T7T and does not conflict with symmetry. Although authors of ref. p0[ | also pointed out that the present experimental 
data still can not be used to separate "direct" uj?:?: coupling from uj — p mixing contribution in model-independent 
way, it is very interesting to perform a theoretical investigation on "direct" ujitit coupling contribution. 

It has been known that uj — > tt + tt~ decay amplitude receive the contribution from two sources: isospin symmetry 
breaking due to u — d quark mass difference and electromagnetic interaction. In section 3 we have shown that VMD 
in meson physics is natural consequence of the present formalism instead of input. The p — > e + e _ decays are also 
predicted successfully. Therefore, the dynamics of electromagnetic interactions of mesons has been well established, 
and the calculation for uj — ► h + tt~ decay from the transition uj — > 7 — * p — > tttt and "direct" uj — ► 7 — * 7T7t is 
straightforward. Then we can determine isospin breaking parameter Sm q = — m u via uj — > 7r + 7r _ decay. This 
parameter is urgently wanted by determination of light quark mass ratios. 



A. Relevant effective vertices at tree level 



Here we first abstract those vertices relating to uj —> irir decay from effective action in section III. 
The p° — uj mixing vertex, which breaks isospin symmetry, is included in eq. (3.13), 



£ ( ) = iV^m^m, 1 _aj_ e - iq . Xa3(q , ){g ^ qtiqv)ul » {q)p o» {x) + 0{{md _ mu fy (6J) 



N c m u — m d f d A q 
6ir 2 g 2 m J (27r) 

The isospin symmetry unbroken vertex-<I><I> vertex (including pnir, AI — pKK and ujKK vertices) and 4-$ vertex 
has been given in eqs. (3.25) and (4.7) respectively. The isospin symmetry broken vertex-<I><I> vertex (including a;7T7r, 
AI = pKK and ujKK vertices) can be abstracted from 5^ normal ) anc j ^| normal ) ^ 

4 A « 1] - aH C J nU ~ md I S^e-^s^ig^q 2 - q,q v )Tr f {X 3 V^q)Hx)d^(x)}, (6.2) 



6tt 2 F 2 m J (2tt 
where A3 = diag{l, —1, 0}, and form factor s(q 2 ) is defined as 

■s(l 2 ) = g-'Mq 2 ) + -/ A (a 7 (q 2 ) - ^)}. (6.3) 

The isospin symmetry broken 4-pseudoscalar vertex is proportional to m c i — m u , which is included in the following 
lagrangian 

= ^fB Tr f {M{U + C/t)} + ^ g 2 A Tr f {V,UV»U\MU^ + UM)}. (6.4) 

It can be found t hat the all relevant vertices are free parameter k independent. Moreover, we can see that every 
vectors ipip in cq.( |6.2| ) includes a n tra nsversal operator (g^q 2 — q^qv) (where g M denotes four- momenta of vector 
mesons). Thus the first term of eq.( 6.4) does not contribute tow-> 7r + 7r~ decay via pseudoscalar meson loops. This 
transversal operator also constrains that the vertices with one of factors KK, d^Kd^K and Kd 2 K do not contribute 
to uj — > 7r + 7r~ decay via pseudoscalar meson loops. Then isospin symmetry breaking KKirir vertices can explicitly 
read as, 

L kk = ^ = ^9 2 A m(m u K+d»K- - m d K°d^K )n + d^- . (6.5) 
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B. Loop correction to "direct" conn vertex 



The loop correction of pseudoscalar mesons to "direct" ujtttt vertex is similar to one to pirir vertex (fig. 2). The 
contribution from tadpole diagram is obtained via integrating quantum fields in w - 4$ vertex, which is abstracted 
from 53 and S4, 



-e lq ' x s{q ){g^q a - g^q v )uj^{q) 
(Tr f {\ 3 (d"n\ a [\ a , d a Tr] + \ a [\ a ,d"Tr]d a Tr)} 



8 

a=4 



+ ^Tr f {I 2 d"nd a ir{\ a \ a M + MX a X a + \ a MX a )} 



= - uXkwk'' 1 "™'' 1 ' 1 I ~ %^^)Tr f {XM^n( X )}, (6.6) 

where I2 = diag{l, 1,0}. 

The "direct" lutttt coupling effective action yielded by AT-loop (fig. 2-b) can be evaluated as 

iS^ lo ° P) = -J ***V < 0\T{C^°\ X )clt I K ^(y) + C$Z%)£££${v)}\° > ■ (6-7) 

The result is 

Cli-J^ - - J C F r U m md J ^e-^6rn^W)[l + -^tf - - s( q 2 )}E K ( q 2 ) 

x (q\u - q^)^{q)Tr f {\^{x)d v Ti(x)}. (6.8) 
The calculation on chain- like approximation correction of pion to lotttt coupling is still similar to one on p — > tvtt decay, 



C. Loop correction to w — p° mixing 

The loop correction of pseudoscalar mesons to uj — p° mixing is similar to one to p-propagator (fig. 5). The 
contribution from tadpole diagram is obtained via integrating quantum pseudoscalar fields in u>p°ipip vertex, which is 
abstracted from cj( normal ) ) 



12n 2 g 2 F 2 m J (2tt) 4 

jwf k 2 - 4 + ie X>/{A 3 (A°A°A< + A4\ a \ a + A°MA°) } 
Xm 2 K N c m u -m d f d 4 q , 2 2 0i/ 



a=4 

e-^a 3 (9 2 )(.9^g 2 -^K(g)p 0iy (x). (6.10) 



6n 2 F 2 6it 2 g 2 m J (2ir) 

The calculations on two-point diagram of if-meson (fig. 5-b) and chain-like approximation of pion meson (fig. 5-c) 
are similar to one in section V. The results is 

4v' oop) - - J c F « ma ~ md I ^ e ^'^^/^>^^ - wy^W, (6.11) 

where T,^\q 2 ) is defined in eq. (fjTg|). 
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D. Electromagnetic correction 



Due to VMD, the electromagnetic interaction contributes tow-* ir + n~ decay through one-photon exchange. The 
"complete" uj — 7 mixing vertex can be obtained via calculating loop contribution in fig. 4-a) and fig. 4-b). It should 
be pointed out that, in this case, the pion loop correction is very small due to suppression of isospin conservation. 
Then we have 



/•(c) 



1 



d 4 q 

W) 4 



-iq-x He) I 2 



m{<t){9^<t - q^)A»(x)Tr f {QV»(q)}, 



(6.12) 



where 



(6.13) 



where T^\q 2 ) is defined in eq. (HH). Then cqs. ( |of ) and ( |6T2|) will lead to uo — p° mixing at order of a c 
through the transition process ui — > 7 — > p°, which is 



/*c.m. 



12 



(27T) 



(6.14) 



Moreover, eqs. (4.31) and (6.12) also lead to "direct" ojtttt coupling at the order of a c . m . through the transition process 
lu — > 7 — > 7r + 7r _ , which is 



/•cm. = J_ 2 

UJTTTT 



d 4 q 



)u,»{q)Tr f {\ 3 Tr(x)d u ir{x)}. 



(6.15) 



E. p° — uj mixing and uj — » tt + tv decay 
Equation ( |6.14 ) together with eqs. ( 6.10| ) and ( |6.11 ) give the "complete" uj — rho° mixing vertex as follow: 



r(c) _ 



(2tt)4 



®u> P (q Z ){ 9 ^q Z ~ q^)uj^q)p 0v (x), 



where 



e wp (g 2 ) = 



N c m u - m d 
6Tr 2 g 2 m 



Am 



W<? 2 )(1 - + 2. 9 2 g 2 /W 7r (< z 2 ) s (g 2 )[s i ,(g 2 ) 



q 2 s£V) 

fIi + e^ 2 



:]} 



The complete "direct" W7T7T vertex can be obtained via summing eqs. Q6.6D , (6.8), ( |6.9D and ( |6.15| ), 



2 7 (2tt) 4 

where the form factor / W7rw (g 2 ) is defined as follows: 



e- tq - x f^(q 2 )(q 2 5^ - q^q v )uj>"{q)Trf{X^{x)d v T:{x)} 



f UJTTTT {.Q ) 



N c m u - m d 
3n 2 F 2 m 



Mi 2 



l + ^(q 2 ) 12tt 2 F 2 



6- 2 4°l(, 2 )[l + ^(9 2 ~ - S(q 2 )) Z K (q 2 )} + ^ < ' 



(6.16) 



(6.17) 



(6.18) 



F^(i 2 )Uy(q 2 )- (6-19) 



Thus G-parity forbidden uj — > tt + 7t inclu des a nonresonant background contribution, cq.( 6.18), and p resonance 
exchange contribution(eqs.( 4.17) and (|6.16|)). The decay width on u mass-shell is 
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487T mz, — m~ + vm^Y p (m^) 

For obtained the mass parameter of p- meson, m p , it is needed to study the difference between mass parameters of 
p° and ui. It is expected the masses of vector meson octet is degenerated at chiral limit and large N c limit. In general, 
the difference between mass parameters of p° and ui is caused by three sources: p° — uj mixing, VMD and loop effects 
of pseudoscalar meson. The loop correction to two-point vertex of p has been obtained in eq. (5.7). Similarly, since 
pion loop is suppressed by isospin conservation, the loop correction to two-point vertex of uj is 

C l J°° P = ~\f ^e- i * x [f^{ q 2 )fml^\ q 2 ) q \g liv q 2 - q^)^(q)^(x). (6.21) 

Furthermore, due to eqs. (4.22), ( |6.12j ) and (6.16), the calculation on VMD correction and p° — to mixing correction 
is straightforward. In addition, since T w <C m^, the mass parameter of uj is almost just the physical mass of to. In 
the other words, m u = m w is a good approximation. Then the final results are 

2 2 , o r ^ ®^p(^ ) , ^.m. 2tr(c)/ 2\i2 4l -c(O) I 2m2 2v-i(0)/ 2\i 

m uj = m v + Re-f— — ^ + — — q [f$(q )f - q [f^(q )) m K T, K K '{q )}| 9 2 =m 2 , 

q* — mr p + q *L p {q z ) ■> 

ml = ml + Re{ ^f^=- + na c . m .q 2 [f^ (q 2 )} 2 
q 2 - m£ + iWq z T u 



vei(? 2 )] 2 ( -K^m (6.22) 



Using the experimental data, m w = 782MeV and B{uj -> tt+tt") = (2.21 ± 0.30)% [ff], together with cq.( |6.20| ), we 
obtain 

m v = 782.5MeV, m p = 803MeV (6.23) 

and 

m u -m d = -(3.87 ±0.21)MeV (6.24) 

at energy scale ~ m u . Here the error bar is from the uncertainty in branch ratio of the process uj — > 7r + 7r~. In the 
standard way, the u> — p° mixing amplitude is 

/ (§F n " p(<z2) =< W| / ^ xC «p( x )\p > => n -p(« 2 ) = 9 2 e^(g 2 ). (6.25) 

The off-shell uj — p° mixing amplitude is obviously momentum dependent, and vanished at q 2 =0. This is consistent 
with the argument by O'Connell et. al. in ref. p4| that this mixing amplitude must vanish at the transition from 
time-like t o spa ce-like four momentum within a broad class of models. In addition, the value of isospin broken 
parameter( 6.24 ) leads on uj mass-shell uj — p° mixing amplitude as follow 

Ren wp (m 2 J ) = -(3888 ± 253)MeV 2 , lmR up {ml) = -(2018 ± 122)MeV 2 . (6.26) 

In ref. pS^ , the on-shcll mixing amplitude has extracted from the e + e _ — ► 7r + 7r _ experimental data in a model- 
dependent way. In eq.( 6.25), the real part of on-shcll mixing amplitude agree with result of ref. ]29| ]. The imaginary 
part, however, is much larger than one in ref. |^9| which is around — 300MeV 2 . It must be pointed out that, in ref. 
p9[ the author's analysis bases on a model without "direct" ujittt coupling. Therefore, it is insignificant to compare 
the value of on-shell mixing amplitude of this the present paper with one of ref. ]29(| . Fortunately, the ratio between 
uj — > 7T7T decay amplitude and p — » 7T7t decay amplitude should be model-independent. This value can test whether a 
model is right or not. The on-shell mixing amplitude in ref. |2S|] yields 

= < n+n J" > = -(0.0060 ± 0.0009) + (0.0322 ± 0.0050)i. (6.27) 

1 < 7T+7T \p > 

The present paper predicts 
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(— ) (tyi 1 f ( in 

R^p = —5 — "I" " M " h 9N - = -(0.0084 ± 0.0007) + (0.0323 ± 0.0021)*. (6.28) 
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We can see that this theoretical pr edict ion agree with experimental excellently. 

If we take /^^(q 2 ) = in eq.( |6 ■ 20[ ) , we have B(ui — ► 7r + 7r~) = (2.16 ± 0.19)%. So that the contribution from 
interference between "direct" lwrtt coupling and u> — p° mixing is very small. The dominant contribution is from 
p-resonance exchange. This result indicates all previous studies which without "direct" cvmr coupling are good 
approximation even though the neglecting the "direct" coupling is an ad hoc assumption. 



VII. to -> 77T° AND p ± -> 7 7r ± DECAYS 



At th e lea ding order of N c 1 expansion, the p'yir and ui^tt vertices are including in anomalous three-point effective 
action ( 3. 16 ) . In momentum space, it is 



n 



N c 



V-yn 



3tt 2 5 F w 



eg A a 6 (q 2 , O)e^q ot k {pUq)A v (k)7r i (q - k) + 3cji(q)A v (k)w°(q - k)}. 



(7.1) 



The loop correction on these anomalous decays are similar to one on p — > kit decay. The contribution from tadpole 
diagram is 



Nr. 



tePgF* 9 * J (2tt) 4 p 2 - m 2 K + it 
x (pi{k)X j + 3uj v (k) + eQA v {k)){[X a ,X l ]X b + X a [X l ,X b }) 
+e[(pl(q)X l + 3w M ( 9 )), ([X a ,Q]X b + X a [Q, X h ])]A v (k)X 1 } 

^Tlvyiriq, k). 



d A p 



id' 



ab 



^ a Pq a k^ l Tr f [^{q)X l + 3^(q) + eQAp(q)) 

(i,j,l= 1,2, 3; a, 6 = 4,. 



Up to 0(p ), the 7-3$ vertex can be abstracted from eqs. ( 3.16 ) and ( 3.21 ), 



n 7 _ 3 $(fci, k 2 , fa) 



^^eg A (l - 2±)eV"*f > k lv k 2a k af }A v (k 1 + fa + k^TrfiQ^faWfaWfa)}. 



(7.2) 



(7.3) 



Then using eqs. ( |3.2E| ) and (7.3), the contribution from two-point diagram of if-meson is 

x{pl(9)4(%.(9 -k) + 3cj i Jq)A v (k)w°(q - k)}. 



(7.4) 



In addition, p^ir vertex also receives contribution from chain-like approximation of pion loops (the contribution to 
Lu^irir is suppressed by isospin conservation). The result is 



TT7T — loOJ) 
P77T 



8tT 2 F, 



■^p'MA^k^q-k). 



(7.5) 



Then eq. ( |7.l| ) together with eqs. ( |7.2| ), (7.4) and (7.5) lead to "complete" p^ir and ujjit coupling to the next to 
leading order of iV" 1 expansion at least, 

Tfy^{q, k) = ee^qMfp^iq^p^A^miq - k) + g^{q 2 )uj,{q)A„{k)^{q - k)}, 



/p77r(<? ) 



3ir 2 gF„ 
N 



9 A (1 



Xm 2 K 

-»2 



)a 6 (q 2 ,0) 



8tt 2 F 
N 



r«.(l-^e(rt( T |^-W)), 



A ^) a6 ( 9 ^o) + ^ 5A (l-^),^/-: r ( ^ ^)V /i( ^ i . 



t 

3 



(7.6) 
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For g A (p = nip) =0.75, the above results yield 

B(p ± -> 71^7) = 4.75 x 10~ 4 , B(u -> vr°7) = 8.8%. (7.7) 

These results agree with data @ B(p ± -> 71^7) = (4.5 ±0.5) x 10~ 4 and B(w -> 7r 7 ) = (8.5 ±0.5)% very well. Since 
in these anomalous decay, the effects of axial constant g A is the leading order, it can be checked whether g A = 0.75 
or not, which is fitted by /3 decay of neutron. If we take g A = 1 as a comparison, the theoretical predictions are 

B(p ± -► 71^7) = 8.51 x 10~ 4 , B{uj -► vr°7) = 15.6%. (7.8) 

These results obviously disagree with data. This point implies that the EMG model are available not only on light 
flavour meson physics but also on light flavour baryon physics. 



VIII. PION FORM FACTOR 



The process e + e~ — > 7r + 7r~ at energies lower than the chiral symmetry spontaneously breaking scale contains very 
important information on low energy hadron dynamics. It was an active subject and studied continually during past 
fifty years. Experimentally, the effects of the strong interaction in process of e + e~ annihilation is obvious to provides 
a large enhancement to production of pions in vector meson resonance region fl37|- 4l| . Theoretically, however, the 
problem was not studied completely so far. In some recent refrences p3|-^5|, the authors have studied e + e~ — > 7r + 7r~ 
cross section and 1 = 1, I — 1 phase shift at vector meson resonance region by using some simple phenomenological 
models. Each of them capture some leading order effects of the chiral expansion and N~ x expansion, and they are 
classified by different symmetry realization for vector meson fields p6| . However, as we shown in previous sections, it 
is not enough to consider the leading order effects only at vector meson scale. Otherwise many delicate features on 
vector meson p hysic s will b e l ost. 

Eqs. J4 1 7D , fl4.22D , fl4.3l| ), (|t|), fl6.16| ) and fl6.18| ) lead to the electromagnetic form factor of pion as follows: 



-^77T7t(Q ) Fry ft -ft (C[ ) 



gV£V)/$r(? 2 ) 



A Ac) / 2\ (c) 

y Ju7w /yum 



(<1 2 



m 2 



(q 2 )) 6(q 2 -m 2 +iy 



(8.1) 



where 



(c) r„2 



(q 2 ) = 



q 2 e^ P f P Uq 2 ) 



_ /(c) 



m + 1* 



(q 2 



f [c> (a 2 ) 



(8.2) 



Here due to narrow width of u>, we ignore the momentum-dependence of T w . In this form factor, we can see that the 
contributions of resonance exchange accompany q factor. Due to this reason, some authors declared that the pion 
form factor in WCCWZ realization for vector meson resonances exhibits an unphysical high energy behaviour(/i > m p ). 
However, this conclusion is wrong. It is caused by their wrong result for momentum-dependence of T p (q 2 ) which is 



fitted by experimental instead of by dynamical prediction. In fact, since 



^p(q ) is proportional to q° 



at least, 



we do not need to worry that the form factor has a bad high energy behaviour. We can also see that there is a 
moment-dependent non-resonant contribution. It together with the contribution of resonance exchange determined 



the high energy behaviour of the factor. The cross-section for e + e 



is given by(neglecting the electron mass) 



(q 2 - Ami) 3 / 2 



\FAq 2 W 



(8.3) 
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FIG. 7. e + e — * -k + tt cross section. The experimental data are from refs. [31,32]. 



From definition of function ai(q 2 ) in section III we can see this model is unitary only for q 2 < 4m 2 . Thus the 
effective prediction should be below < 2m = 960MeV. The result is shown in fig. 7. We can see the prediction 
agree with data well. Especially, the theoretical prediction in vector meson energy region agree with data excellently. 
Although the mass parameter rh p = 803MeV in p propagator is larger than physical mass, the position of pole is 

localized in y/q 2 — 772MeV which is just the physical mass of p. It strongly supports our above discussion and 
dynamical calculation. It also implies that we must carefully distinguish the physical mass difference of p° and u> from 
the difference of mass parameter in effective lagrangian. 



Let us give some further remarks on pion form factor (3.2). In eq. (S.2), the form factors F 1 ^ 1T {q 2 ), f pl (q 2 ), etc., 



are all complex function instead of real function. It is caused by one- loop effects of pions. Thus the expression (8.2) 
can be rewritten as follows: 

2{q z -171*+ ly/q^L p (<r)) b{q z — m£ + i y q 1 u) 

Here bi(q 2 )(i = 1,2,3) are three real function and </>i(g 2 )(« = 1,2,3) are three momentum-dependent phases. In 
particular, </>3 (q 2 = m 2 ), so called Orsay phase, has been extracted from data as 100 — 125 degrees Hj45||. Our 
theoretical prediction is 4>?,{q 2 = m^) = 116.5 degrees. However, so far, the phases 4>\{q 2 ) and 02 (tp) are not 
reported in any literature. These momentum-dependent phases indicate that the dynamics including loop effects of 
pseudoscalar mesons is different from one only in tree level. In fig. 8, we given theoretical curves of 4>i(q 2 ). They are 
indeed nontrivial. 
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FIG. 8. (pi versus m-mr in GeV. Here the solid line denotes the phase shift of non-resonant 
background tj>i, the dash line denotes the phase shift 02 in p coupling and the dot line 
denotes the phase shift 03 in uj coupling. "•" denotes Orsay phase. 



Obviously, ^(q 2 ) is an analytic function in the complex q 2 plane, with a branch cut along the real axis beginning 
at the two-pion threshold, q 2 = Am 2 . Time- reversal invariance and the unitarity of the S- matrix requires that the 
phase of the form factor be that of I = 1, I = 1 ir — it scattering p7| . This last emerges as ir — ir scattering in the 
relevant channel is very nearly elastic from threshold through q 2 ~ ( m^ + m u ) 2 Q^^]. In this region of q 2 , then, the 
form factor is related to the 1 = 1, I = 1 w — tt phase shift, S\, via EjJ 

F^{q 2 )=e 2 ^F;^{q 2 ), (8.5) 

so that 

tan5l(g } = ReF^Y (8 - 6) 

The above is a special case of what is sometimes called the Fermi- Watson- Aidzu phase theorem jl9|,^0|]. In fig. 9 
and fig. 10 we plot theoretical curves of the 1 = 1, I = 1 tt — ir phase shift 5{ versus m W7r and of sin 5{/p^ versus 
(where = \\/ q 2 — Am 2 ) respectively. We omit the u> contribution from our plots of the phase of F n (q 2 ) for 
comparing with time-like region pion form factor data ]5l|-|5l| . We have also assumed that S\ is purely elastic in the 
regime shown, i.e., the loop effects of u> — ir are omitted. The curve predicts <5j — > 90° as y^q 2 — » 774MeV ~ m p , and 
8\ > 100° for y/q 2 > 787McV. These results agree with data very well. 
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FIG. 9. The 1 = 1, I = 1 rr — n scattering phase shift. The solid circle point are these 
from [44], the hollow diamond point are from [45] and "+" denotes the points from [46]. 
Note that the uj contribution to the time-like pion form factor phase has been omitted, to 
facilitate comparison with the empirical phase shift. 



Finally we discuss the near threshold behaviour of the form factor. 1) The chiral perturbative theory predicts the 
form factor at threshold to be [F 77rw (4m^)]chPT = 1.17 ± 0.01, and ours, [F 7r (4m^.)] = 1.145, is close to the ChPT 
result. 2) The electromagnetic radius of charged pion has been determined to be \f< r 2 > v = 0.657 ± 0.027fm Q, 
whereas the theoretical prediction in this present paper is \f< r 2 > n = 0.645 fm. 3) The Froggatt-Petersen phase 
shift function sin 5{/p^ is connected with the vector-isovector it — it scattering length a\ through 

i , sin^, 1 

a\ = lim — (8.7) 

Our theoretical prediction is a\ = 0.037 in unit of m~ 3 . This value is very close to experimental results from K e 4 
data fUJiH usin S a R °y equation fit (a{ = 0.038 ±0.002) and ChPT prediction a\ = 0.037 ±0.01 (57) at the two loop 
order (at 0(p 4 )). 

To provide a brief summary on this section. In section V, it has been revealed that the mass parameter in resonant 
propagators should be different from its physical mass due to large momentum-dependent width. In this section, this 
point is confirmed by both of dynamical calculation and phenomenological fit. It also tells us how to understand 
the mass splitting between p° and to: Although the dynamical calculations show that the mass parameter of p in 
effective lagrangian is even larger than one of w, the position of pole localized in real axis give their right physical 
mass splitting. The contribution to this mass splitting from p° — u> mixing is very small, the dominant contribution is 
from one-loops effects of pseudoscalar mesons. The effective field theory mechanics on p° — u) mass splitting revealed 
in this present paper is rather subtle. And it is another evident to confirm again that the EMG model is sound. 
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IX. LIGHT QUARK MASSES 



The quark masses are some of basic parameters of the standard model. There are various hadronic phenomenology 
relating to the light quark(u, d, s) masses. For instance, they break the chiral symmetry of QCD explicitly. m u — md 
breaks isospin symmetry or charge symmetry, and (m u + m^jl — m s breaks SU(3) symmetry in hadron physics 
respectively. However, in QCD the masses of the light quarks are not directly measurable in inertial experiments, but 
enter the theory only indirectly as parameters in the fundamental lagrangian. Therefore, it was an active subject to 
determine light quark masses via phenomenological method. 

At low energy, the information about the light quark mass ratios can be extracted in order by order by a rigorous, 
semiphenomenological method, chiral perturbative theory (ChPT) To the first order, the results, m s /m c i = 19 

and m u /md — 0.556, are well-known fl59"| . Many authors have studied the mass ratios up to next to leading order 
of the chiral expansion. Gasser and Leutwyler first obtained m s /md — 20.2 and m u /md — 0.554 Q. Then Kaplan 
and Manohar extracted m s /md =15 to 23 and m u /md =0 to 0.8 with very larger error bar |60|. These values have 
been improved to m s /md = (20.5 ± 2.5), m u /md = 0.52 ± 0.13 by Leutwyler [plfl, to m s /nid = 18, m u /irid = 0.66 
by Gerard p^ , and to m s /md = 21, m u /md = 0.30 ± 0.07 by Donoghue et al. ]63| respectively. Finally, Leutwyler 
analysis previous results and obtained m s /md — 18.9 ± 0.8, m u /md = 0.553 ± 0.043 fl64f . 

Although ChPT provides many information on light quark masses, it is still not enough for further studies on 
vector meson physics: At chiral limit, the unitarity implies that EMG model can not predict the physics when energy 
yfip > 2m ~ 980MeV. In the other words, to study ^(1020) physics at chiral limit is not consistent. However, if we 
take nonzero light quark masses, EMG model will be unitarity when energy > 2(m + m s ). It means that all light 
flavour vector meson resonances can be included in to EMG model consistently. Hence, the strange quark mass plays 
important role at K* and <j) physics. Meanwhile, ChPT can only provide information on light quark mass ratio. For 
obtaining the individual quark masses, other approaches, such as QCD sum rules |65||36| or lattice calculation |67| are 
needed. Some authors have pointed out that the determination on the individual quark masses is model-dependent 
|p8| . Thus for studying K* and <\> physics in formalism of EMG model, it is better to determine light quark masses 
by this model itself. 
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The Kaplan-Manohar ambiguity of ChPT has caused many debates. In particular, due to this ambiguity, the 
authors of ref. |(38f| argued that the observed mass spectrum is consistent with a broad range of quark mass ratios, 
which specially includes the possibility m u = 0. Howe ver, it is disagreed by other anthors p9,p3[. This problem will 



also be discussed in EMG model. From lagrangian (2.12), we can see that light current quark masses are defined 



uniquely in this model, that they are just renormalizcd "physical" masses of u, d, and s quarks. In principle, therefore, 
there is no Kaplan-Manohar ambiguity in EMG model. Light current quark masses can be determined uniquely via 
meson spectrum. 

In general, in EMG model the light quark masses can be extracted not only by pseudoscalar meson spectrum, 
but also by lowest vector meson resonance spectrum. However, one-loop effects of mesons will contribute to vector 
meson masses. As shown in the above sections, the calculations to every vector meson resonance are very complicate. 
It makes that the relationship between vector meson spectrum and the light quark masses are also very complicate 
and indirect. Thus here we will extract information about the light current quark masses from pseudoscalar meson 
spectrum and their decay constants. The vector meso n spectrum will be predicted by this formalism in other paper. 

For the purpose of this paper, the lagrangian ( [2.12 ) can be rewritten as follow 



C x (x) = C q (x) + 4 0) (z) + C[w(x),r]s(x)} + £[if ± (x)] + C[K°(x)], (9.1) 
where C q and are f ree field lagrangian of constituent quarks of pseudoscalar meson respectively, 

Cq = ^ qi(ifl - fhi)qi, fhi = m + mi, 

F 2 

= —d fl <S> a d> 1 <S> a + 4a a ^<S> a + ---, a = 1, 2, • • • , 8. (9.2) 



jC[$(x)] denotes quark-meson coupling lagrangian, 



C[n, m ] = %(d M 7r+ + 2a+)u 7 " 75 d + «.] + 9 -f ]T (d^P + 2Af )q t 



7^75 1i 



i—u,d,s 



H — j=K{m u + md)(7r + uj^d + c.c.) +itz ^ m i P l q i j^q t 

i—u,d,s 

+ 2^ mu + m d )7T + 7T~(uu -t- dd) + i ^ TmFfqi9u (9-3) 

i=u,d,s 

£[#+] = ^[(d^K+ + 2A+)u^ l5 s + c.c] + A^K(mu + m s ){K + u l5 s + c.c.) 

+ 7j( m « + m s )K + K~(uu + ss), 

C[K°] = ^[(d,K° + 2A°)dy 75S + (d,K° + 2Al)d^ 75 s] 

i - - 1 

+ —=n{m d + m s )(K°dj 5 s + K°dj 5 s) + -{m d + m s )K° K a (dd + ss), 
V 2 2 

where c.c. denotes charge conjugate term of previous term, 

1 12 , A 

P u = tt 3 + —^=r] S , Pd = -iT3 + -^m, ^^""Ta 778 ' ^ 9 ' 4 ^ 

and o°, A^\ A/?' and A^ etc. are axial- vector external fields corresponding these pseudoscalar meson fields. From 



eq. (EL3) we can see that the -ft^-quark coupling and i4f°-quark coupling are similar to ^-quark coupling. Thus we 
only need to calculate masses and decay constants of , 7r° and rf . Then masses and decay constants of K ± can be 
obtained via replacing by m s in one of tt^ , and masses and decay constants of K° can be obtained via replacing 
m u by m s in one of . 

At large N c limit, the kinetic terms of pseudoscalar mesons can be obtained via calculating one-point and two- 
point effective action. Note that here the light quark masses will not be expanded perturbatively, but enter quark 
propagator and be evaluated in non-perturbative method. Explicitly, the kinetic term of charge pion fields reads 



2cS 



where 



3 

F 2 (m u , m d ) = Fq + -g 2 g 2 A m u d(2m + m ud ) + ig 2 ng A m u drh ud 

--^9A(gAm u d + 2Km u d) J dt[m u + t(m d - m u )] In {x 2 u + t(x d - x 2 u )) 

, , 2 2 U 9 2 F B N c s N c f 1 2 2 2 m„m d n 

+4k ^{(—--—,-—5)-—^ / dt-t(l-t)[31n(g u +t(g d -gj)- - , A . 2 — ^2T]}' 



4 48m 3 247T 2 ' 8tt 2 ,/„ v n v u v <* m^^.^)! 

rtm.,md - F„* + §.««*«-(«» + *-) + ^W-ffw 

~~2~2 gA ( 9Af?lud + K17lud + y *K + t ( m d ~ m «)l ln ( S m + - ^u))' 

M 2 (m u ,m d ) = B m ud {hxl + x%) + ^ c ™ (x 3 u \nxl + x d \nx 3 d ) - — (x 2 + a; 2 ,)} 

3 g2/t 2 K ~<) 2 _ .^^^L / dt[ 2rhl-rh u rh d + 2t(m 2 -rhl)}\n(x 2 u + t(x 2 ~x 2 u )) : (!).(>) 

,2. 



2" 2 7T 2 F Q 2 Jo 



iV c , , Z" 1 , r _ , .„ / t(\-t)q 2 
a(q ;m u ,mdj ^ -^^gA(gAm u d + Km u d) / at[m u + t(m<i - m„)J ln ( 1 - 



2tH 7 V TO d + H m d- m «)/ 

i«, md) = 7T^9A-^2(gAm u d + 2K,m ud ) / di[m u + Z(m rf - m„)] In 1 - _ 2 . _ 2 r^r 

z!7r ^ Jo \ m u + t(m d — m u ) 

^n 2 ^$f t dt{[2fn 2 u - m u m d + 2t{m 2 - m 2 ) - 3t(l - t)g 2 ] In f 1 - ry j^^ 
- m~ + t{m d — to„) 



•0 ^0 



2- . ^ }, 
V ™„ + H™d-™u)/ 

where my = nij + mj and my = m, + fhj(i,j — u,d,s). It should be pointed out that a(q 2 ;m u ,m d ) is order q 2 
at least and /3(q 2 ;m u ,rrid) is order g 4 at least. Since in this paper we focus on pseudoscalar meson spectrum, these 
high order derivative terms should obey motion equation of pseudoscalar mesons. In momentum space, the motion 
equation of physical pseudoscalar mesons is generally written 

(q 2 - mlMq) = -if^ A^\-q), (9.7) 

where m v and f v are physical mass and decay constants of pseudoscalar, e.g., = 135MeV and f n — 185.2MeV. 
Due to this motion equation, we have 

a{q 2 :m u ,m d )a+(x)ir~(q) = a(m 2 ; m u , m d )a+(x)ir~ (q), 



f3(q 2 ;m u ,m d )Tr + (q)TT (-q) = /3(m 2 ; m u , m d )n + (q)TT (-q) - ^^/^^'(m 2 ; m u , m d )(a+(x)n (q)+c.c), (9.8) 



where 

(m 2 : m. . . m. j \ = 

dq 



/3'(m 2 ;m u ,m d ) = — f3'(q 2 ; m u , m d )\ q 2 =m 2 . (9.9) 



Thus eq.( pM can be rewritten 



£ 2 [tt±(x)} = F2{m 2 ,md) d^+d^- - ^{M 2 (m M , m d ) + /3(m 2 ; m u , m d )}ir+ir- 

r f 2 (m„,,m d ) a(m 2 \ra n ,m d ) ,,9 _i_„„ 

+{ V f ' + 1 - 2 ' + 4F{m ° md) UP'{ml- m u ,m d )}{a+d^- + c.c). (9.10) 
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The kinetic term of 7r° and r/g reads 
where 



F g F Bq 

F?{mi) = + 3g 2 g 2 A m t (2m + m l ) + 6g 2 ng A m i m l + 8k 2 (— - , 
2 4 48m J 48ir z 



e i «- x /3 i (g 2 ;m i )Pi(g)Pi( i r)} ) (9.11) 



A, 



•^c / 2 - 2 , n - , 2 2\ i 2 

- ^2 (ff A m » + ^gAKUiirrii + n m J In ^ , 



M 2 (mi) = (S 
ft(? 2 ;m.) 



iV c TO 3 TOj 

— -p2 hiiTjJmiZj^i - — k ), 

7T _T n 



m 



(9.12) 



2F2 g 2 {^- m » + (gA^i + Znrrufhi + 2K 2 TO 2 [m 2 g~ 2 - 3t(l - t)]) / diln(l- 
o ^ Jo 



£(1 — t)g 2 



)}■ 



Since the decay constants for neutral mesons cannot be extracted directly from the data. It means that the decay 
constants for neutral mesons can not be used to determine light current quark masses. Therefore, in this paper we do 
not need to evaluate the decay constants for neutral pion and rjg. In addition, because auxiliary fields Pi do not lie 
in physical hadron spectrum, the equation of motion (9.7) ca n not be used in (9.11) simply. We will use propagator 
method to deal with the terms with high power momenta in (9.11) and diagonalize — rjs mixing. 

Next we evaluate one-loop effects of pseudoscalar mesons. Due to parity conservation, there are only tadpole 
diagrams of pseudoscalar mesons contributing to masses of decay constants of 0~ mesons(fig. 1). Moreover, we can 
expect that, in mass spectrum and decay constants of pseudoscalar mesons, the dominant one-loop effects is generated 
by the lowest order effective lagrangian, because neglect is 0{m^) and suppressed by A"" 1 expansion. It is fortunate 
that the formalism is renormalized up to this order. 

The lowest order effective lagrangian is well-known 

Fn 



16 

= H < [A a , A^] [A a , A M ] > +|| < A-A^fc^ + ^t) >, 



(9.13) 



Substituting expansion ( 4.1 ) into £2 and retaining terms up to and including ip 2 one obtains 

£2 - £2 + ^(fl^'W - m 2 (^V) - H < [<p, AJb, A"] + <p}(&t + f " 2M) >, (9.14) 

where we have omitted some terms which do not contribute to masses and decay constants via one-loop graphs. The 
contribution of tadpole graphs can be calculated easily 



£ 



(tad) 



-\ < [A Q ,A Al ][A a ,A^] + i{A Q ,A a }(ex t e + e t X? t -2X) > 



(2tt) 4 k 2 - m% + ie 



16tt 2 



In- 



-) < [\ a ,A M ][\ a ,M] + 7 {\ a ,\ a }(&i + ?X?-2M) >, 



(9.15) 



where 



N, 



167T 



■{- + ln (4tt) + 7 + 1}. 



Comparing eq.( |9.15 ) and eq.( 9.12 ) we can see that the divergence N e can be absorbed by free parameters F and 
Bq. Thus the sum of tree graphs and tadpole contribution is 



AD 



F 



ZF 2 



4 |(1-3 M > < [A",A / J[A» ) A/'] > +2^(1- 3^) < A a A Q (e x t C + ^t) >; 



(9.16) 



where 
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in 



4 ^2 F 2 111 „2 



Here we appoint that m v = m„ when a = 1,2,3, m ¥ 



when a 



(9.17) 



4, 5, 6, 7 and m„ 



'■')s 



when a = 8. 



For extracting (m + )qcd from experimental data, the electromagnetic mass splitting of if-meson is required. The 



prediction of Dashen theorem |7C 



K + 



1.3MeV, has been corrected in several recent analysis with 



considering contribution from vector meson exchange. A larger correction is first obtained by Donoghue, Holstein 



and Wyler |71[, who find (m 



2.3MeV. Then Bijnens and Prades 173] , who evaluated both long- 



distance contribution using ENJL model and short-distance contribution using perturbative QCD and factorization, 



find (m K+ -m K0 ) e . r 



2.4±0.3MeV at 



Gao efc.ai. J73| also gave (m 



l K o )e.r 



2.5MeV. Baur and Urech 



however, obtained a smaller correction, (ni K+ — m K0 )e.m. = 1.6MeV at /j = m p . In addition, calculation of lattice 



QCD J7J found (m jfH 



1.9MeV. These estimates indicate that the corrections to Dashen theorem arc 



indeed substantial. In this the present paper, we average the above results and take (m _, 
at energy scale fj, = m p . 



)e.r. 



2.1±0.1MeV 



From eqs. (9.10) and (9.16), the masses and decay constants of koan and charge pion can be obtained via solve the 



following equations 



in 



in 



/tt-i 



' K+ 



f K 





F 2 (m u , 


m d ) 


F Q 2 




Fl{m u 


m s ) 


F 2 




F 2 (m d , 


m s ) 




™d) 


F R (m u ,m d ) 


fl{m u , 


m s ) 


F R (m u 


m s ) 


Il{m d , 


m s ) 


F B (m d , 


m s ) 



{M 2 R {m u ,m d ) + (3(ml;m u ,m d )} 

{M 2 R {m Ul m s ) + l3{m 2 K ;m u ,m s )}, 

{M R (m d ,m s ) + P(m 2 K ;m dl m s )}, 
a(ml;m Ul m d ) 



(9.18) 



Fi 



F R (m u ,m d ) 
a{m 2 K ;m Ul m s ) 

F R (m u ,m s ) 
a{m 2 K ;m dl m s ) 

F R (m d ,m s ) 2F 2 (m d ,m s ) 



2F 2 (m u ,m d ) 
F 2 

2F 2 (m Ul m s ) 



f w +P'(ml;m u ,md) 



Fi 



f K P'( m2 K ; m u,m s ), 



where subscript "R" denotes renormalized quantity, 






F 2 (m u ,m d ) 


= F 2 (m u ,m d ) - Fq(2hk 






F 2 (rrii,m s ) 


3 

= F 2 (mi,m s ) - -Fg 2 ^ 


f 2fi K +fi ris ), 




M R {m u ,m d ) 


= M 2 (m u ,m d ) - B (m u 


^m d ){-^ + n K 4 




M 2 R ( mi ,m s ) 


= M 2 (m il m s ) - B Q (mi -+ 


3 3 
-m s ){-^ + -n K - 


5 


P R (mu,m d ) 


= f 2 (m u ,m d ) - Fl(2^ - 






f 2 R {m u m s ) 




Y2[i K +(J, Vs ), 





(i = u, d) 



(9.19) 



Here the quantity /x v depends on the renormalization scale fi. It has been recognized that the scale plays an important 
role in low-energy QCD. In this formalism, many parameters, such as light current quark masses, constituent quark 
mass m, axial- vector coupling constant g^, are scale-dependent. The characteristic scale of the model is described by 
the universal coupling constant g which is determined by the first KSRF sum rule at energy scale /i = m p . Hence we 



take fi 



770MeV in pb v , and the physical masses and decay constants, however, should be independent of the 



renormalization scale. 

In ord er to obtain the masses of tt° and rf , the ^3 — 778 mixing in eq. (9.11) must be diagonalized. Eq.( |9.11 ) together 
with eq.( |9.16| ) lead to the quadratic terms for the ir 3 and rj 8 are of the form 



S^zM = J ^i\(l 2 M 2 3 ( q 2 ))n 2 3 + \{q 2 M 2 8 (q 2 )) V 2 8 - M 2 38 (q 2 )n 3 n 8 } 



(9.20) 



where 
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p2 3 1 

Ml(q 2 ) = -j^{M 2 a {m u ) + M 2 d {m d ) + f3 u (q 2 ; m u ) + l3 d (q 2 ;m d ) - 2B m(-^ + M K + g/^J}. 

M sV) = ^L{M 2 u {m u ) + M 2 {m d ) + AM 2 {m s ) + f3 u (q 2 ;m u ) + p d (q 2 ;m d ) + Af3 s {q 2 -m s ) 

2 3 1 

-2B (rh + 2m s )(2n K + -fi m ) - 2B m(--n„ + n K + -fi Vs )}, (9.21) 

F 2 - - 3 1 

Mls{q 2 ) = - {M 2 (m u ) - M 2 d {m d ) + p u (q 2 ;m u ) - (3 d (q 2 ;m d ) - B (m u - m d )(-^ + n K + -^ 8 )} 



2 F 2 (m u ) - F 2 (m d ) 
V3F 3 F 8 

with m = (m u + m d )/2 and 

F 2 = F 2 (m u ) + F 2 (m d ) - F 2 (2^ + M J, 

*8 = \{Ft(mu) + F 2 (m d ) + AF 2 {m s )} - 3F V- (9-22) 

Due to 7T3 — 7/8 mixing, the "physical" propagators of tt° and rjs are obtained via the chain approximation in 
momentum space 

i * , ^W) .... 



-?2 7T7 ^ 



m^+ie q 2 -M 2 {q 2 )+ie (q 2 - M 2 {q 2 ) + ie)(q 2 - M 2 (q 2 ) + ie) 2 

i 



q 2 -m 2 +ie q 2 - M 2 {q 2 ) + ie (q 2 - M 2 (q 2 ) + ie)(q 2 - M 2 (q 2 ) + ie) 2 



+ ■ 



n 2 _ M 2 (n 2 ) - M 3*( q2) +j f ' 
1 1VI 3 W ) 9 2_ M i (9 2) + «e 

Then the masses of tt° and r/s are just solutions of the following equations, 



(9.23) 



mio = M 3 2 (r<o) 



ml a -M 2 {ml o y 



K- = ^K,) + _ M| '» m (9-24) 



In the above results, the parameters n, Fq and /3o are still not determined. In order to determine them and three 
light quark masses, six inputs are required. In this paper we choose f n + = 185.2 ± 0.5MeV, / + = 226.0 ± 2.5MeV, 
m^o = 134.98MeV, to k0 = 497.67MeV and (to k+ ) QC d = 491.6±0.1MeV. Another input is m d -m u = 3.9±0.22MeV, 
which is extracted from u> — ► ir + ir~ decay at energy scale fi = m p . Recalling m — 480MeV, qa = 0.75 and g = tt^ 1 
for iV c = 3, we can fit light quark masses as in tabic 1. 

In table 1 , the errors of results are from uncertainties in decay constant of K + . electromagnetic mass splitting of 
if-mesons and isospin violation parameter m d — m u respectively. The first column in table 1 we show the results for 
/ + = 223.5MeV. The third column corresponds the center value of / + , 226.0MeV and the fifth column corresponds 
f i<+ = 228.5MeV. Then from table 1 wc have 

m s = 160 ± 15.5MeV, m d = 7.9 ± 2.7MeV, m u = 4.1 ± 1.5MeV, 

fry ryi 

— = 20.2 ±3.0, — = 0.5 ± 0.09. (9.25) 

m d m d 

Here the large errors are from the uncertainty of / + . We can see that our results agree with one obtained by other 
various approaches(e.g., QCD sum rule or ChPT) well. From table 1 we also have (m T + — m w o)QCD = — 0.25MeV 
(where the contribution from tt^ — rj' mixing is neglected) . This result allows that the electromagnetic mass splitting 
of pion is {m^+ — m 7r o) e . m . = 4.8MeV. 
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Finally, it is interesting to expand non-perturtation results (5.1) up to the next to leading order of the chiral 
expansion and compare with one of ChPT. If we neglect the mass difference m d — m u , up to this order the masses of 
decay constants of koan and pion read 



,2 



uf , ui , 1 1 , r 1 ^ 2n , 3m 2 , m g A B 2 m u +m d 
B {m u + m d ){l + ~n* - -fJ, m + I- (3 - k ) ± (— -Kg A - — - —g A )\ }, 

Z Z 7T _Tq Jjq * 0772 771 

,2 ™ „2 



9 „ , « Nr 1 rl/ 9x 3777 ,771 q A Bn n „ , 777 + 777 , . 

m 2 . = Bo(m + m a ){l + + [-(3 - K 2 ) + " ~ f - ^Sa)]— ( 9 - 26 ) 



» c, ri 1 3 2 m(m„ +m d ) 

/tt = ^b{l - Mtt - 2^*- + ^2 h 



3 3 3 3 2 777(777 + 777 s ) 

Comparing the above equation with one of ChPT [Q, we predict the 0(p 4 ) chiral coupling constants, £4, L5 Lg and 
Lg, as follow 

3m 2 



£4 - ie - 0, L 5 - 32^5a 



2 J 

■Fn , On 3777 , 777 Q 2 , Bn n „ L5 , 

;(3-« 2 ) + 77-T7r(ir-^-^-^i) + ^. (9.27) 



128B 777 v ; 647r 2 S v S a 2 6m yA; 2 

Numerically, inputting k = 0.35 ±0.15, F = 0.156GeV and B = 1.7±0.4GeV, we obtain L 5 = (1.2 ±0.3) x 10~ 3 and 
L s = (0.0.47±0.25)x 10~ 3 . These value agree with one of ChPT, L 5 = (1.4±0.5) x 10~ 3 and L s = (0.47±0.3) x 10~ 3 
at energy scale fi — m p . This result indicates that the contribution from scalar meson resonance exchange is small. 
In fact, in hadron spectrum, there is no scalar meson octet or singlet which belong to composite fields of qq. Thus it 
is a ad hoc assumption to argue some low energy coupling constants, such as L§ and Ls, receiving large contribution 
from scalar meson exchange. 





Fit 1 


Fit 2 


Fit 3 


Fit 4 


Fit 5 




185.2 


185.2 


185.2 


185.2 


185.2 


J K+ 


223.5 


224.1 


226.0 


227.8 


228.5 




134.98 


134.98 


134.98 


134.98 


134.98 


if 


497.67 


497.67 


497.67 


497.67 


497.67 


( m K+ )QCD 


491.6 ±0.1 


491.6 ±0.1 


491.6 ± 0.1 


491.6 ±0.1 


491.6 ±0.1 


B(ui — » 7T7r) 


1.95% 


1.95% 


(2.11 ±0.20)% 


(2.21 ±0.25)% 


(2.21 ±0.30)% 


m s 


144.2 


150.0 


161.8 


172.0 


175.4 


m d 


6.26 


6.83 


7.94 ±0.07 


8.94 ±0.1 


9.53 ±0.11 


m u 


2.56 


3.13 


4.10 =p 0.07 


5.03 =p 0.1 


5.63 =f 0.11 


m u + rn d 


8.92 


9.96 


12.04 


13.97 


15.16 


m d - m u 


3.7 


3.7 


3.84 ±0.14 


3.9 ±0.2 


3.9 ±0.22 


f K o 


224.0 


224.7 


226.6 


228.4 


229.0 




134.74 


134.74 


134.73 


134.73 


134.73 


m m 


608.8 


616.9 


628.8 


636.8 


641.0 


K 


0.5 


0.4 


0.3 


0.25 


0.2 


F 


156.95 


156.7 


156.1 


155.6 


155.3 


B 


2141.1 


1886.6 


1551.1 


1328.6 


1216.8 



TABLE II. Light current quark masses predicted by the masses and decay constants of 
pseudoscalar mesons at energy scale [i — m p . Here k is dimensionless, other dimensionful 
quantities are in MeV, and "t" denotes input. 
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X. SUMMARY AND DISCUSSION 



This paper is motivated by two purposes: Theoretically, convergence of chiral expansion at vector meson energy 
scale should be studied. Such a study is important, since it determines whether a well-defined chiral theory can 
be constructed in this energy region. Phenomenologically, we hope to study vector meson physics in more rigorous 
way. In the other words, we hope that these studies are not limited at the leading order of chiral expansion. For 
achieving these two purposes, we need first choose a model, which should satisfy the following requirements: 1) This 
model must yield convergent chiral expansion at vector meson energy scale. 2) Unitarity of S'-matrix in this model 
must be insured. 3) This model must be available to evaluate high order contributions of low energy expansion. 4) 
The theoretical predictions can be derived practically, and the results match with data. Obviously, the ChPT-likc 
models p6| , [76| , [77| do not satisfy the third requirement, since too many free parameters are needed as perturbative 
order raising. In addition, the Extend NJL-like models do not satisfy the first and the second requirements at 
vector meson energy scale. Therefore, EMG model, which satisfy all above requirements, is the simplest and practical 
model for achieving the purposes of this paper. 

There are usually three perturbative expansions in energy below CSSB scale. They are external momentum expan- 
sion, light quark mass expansion and iV" 1 expansion. It is well-known that the light quark mass expansion and N^ 1 
expansion are well-defined for this whole energy region. In the other words, in these two expansion, the dominant 
contribution is from the leading order of perturbative expansion and the calculation for this order is a controlled 
approximation in that error bars can be put in the predictions. This point has been explicitly confirmed in this paper. 
For example, the direct calculation show that high order contribution of N^ 1 expansion is indeed about T p : m p 
which is around 20%. However, the external momentum expansion is not similar to other two expansions. According 
to usual argument on external momentum expansion, this expansion should be in powers of p 2 /A 2 ,. For p 2 ~ m 2 p , 
it is obvious that the momentum expansion converges slowly at vector meson energy scale. In the other words, the 
leading order contribution is not dominant, so that it is very incomplete to study vector meson physics only up to the 
leading order of momentum expansion. In this paper, we show that the high order terms of momentum expansion 
indeed yield very large contribution for vector meson physics. It implies that a well-defined model including vector 
meson resonances must be available to evaluate high order contribution of momentum expansion. 

Some standard methods, such as heat kernel method, are difficult to capture high order contribution of momentum 
expansion. In this paper, therefore, we propose the proper vertex expansion method to deduce effective action. The 
resulting effective action includes all order information of momentum expansion of vector mesons. Hence, in this 
paper, we have achieved a complete study on the chiral expansion theory of EMG model at vector meson energy scale. 

It is crucial and challenging to prove the unitarity of S'-matrix for low energy effective theories of QCD. This 
problem has been solved satisfactorily in this paper. The unitarity has been proved order by order in iV" 1 expansion. 
It yields a strong constrain for all effective theories of low energy QCD, that the transition amplitude must be real at 
the leading order of -/V" 1 expansion. This constrain can be used to examine whether an effective theory is eligible or 
not. 

In this paper, we study a number of phenomenological processes successfully. It should be pointed out that our 
studies arc different from one in all previous literature, since high order contributions of momentum expansion and iV" 1 
expansion arc included in ours consistently. Some new features are revealed in this papers: 1) For wide resonances, the 
mass parameter in its propagator is in general different from its physical mass. However, the difference between them 
is model-dependent. So that the masses of resonances enter effective lagrangian only as a parameters. Their physical 
masses should be predicted by relevant scattering processes. 2) In oj — > n + n~ decay, a nonresonant background is 
allowed, but the contribution from this back ground is very small. 3) There are three nontrivial phases in pion form 
factor (fig. 8). All of these phases are generated by pion loops and are momentum-dependent. 

Up to the next to leading order of A*"" 1 expansion, all results in EMG model are factorized by F n , Bo, m v , g A = (0.75, 
(3 decay of neutron), g = (tt^ 1 , KSRF (I) sum rule), ra=(480MeV, chiral coupling constant at 0(p 4 )), A=(0.54, Zwcig 
rule) and masses of pseudoscalar mesons. So that EMG model provides powerful theoretical predictions on low energy 
meson physics. Furthermore, since strange quark mass has been extracted from pseudoscalar meson spectrum, the 
calculation of this paper can be easily extend to cases including K*(892) and 0(1020). These studies will be found 
elsewhere. 
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